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WAVE  PROPAGATION  ACROSS  ION  CYCLOTRON  RESONANCE  HARMONIC  LATERS 

a ) 
Kaya  Imre    and  Harold  Weitzner 

Coarant  Institute  of  Mathematical  Sciences, 
New  York  University,  New  York,  New  York  10012 

ABSTRACT 
Wave  propagation  across  ion  cyclotron  resonance  layers  is  investigated 
by  use  of  a  boundary  layer  analysis  for  a  perpendicularly  stratified,  weakly 
inhomogeneous  plasma.  Corresponding  to  the  situation  in  many  plasmas  of 
interest,  it  is  assumed  that  the  wave  frequency  is  much  less  than  the  ion 
plasma  frequency,  but  that  the  ratio  of  ion  plasma  frequency  to  wave  fre- 
quency times  the  ratio  of  ion  thermal  speed  to  the  speed  of  light  is  small. 
Non-trivial  effects  of  the  cyclotron  resonance  appear  in  the  boundary  layer 
and  modify  the  current  density  from  the  cold  plasma  limit.  These  correc- 
tions suggest  the  breakdown  of  the  geometrical  optics  limit.  Wave  equations 
for  the  field  valid  within  an  arbitrary  harmonic  resonance  layer  are  de- 
rived, as  well  as  for  the  two  ion  hybrid  resonance  case.  Analytical  proofs 
of  reciprocity  relations  are  given  for  the  general  wave  propagation  prob- 
lems. For  fundamental  resonance,  second  harmonic  resonance,  and  ion-ion  hy- 
brid resonance  proofs  of  the  properties  that  (i)  the  transmission  coeffi- 
cient remains  the  same  as  its  geometrical  optics  value,  and  (ii)  the  absence 
of  reflection  of  the  fast  mode  incident  from  the  high  field  side  are  given. 
Numerical  results  for  the  cases  of  fundamental,  second,  and  third  ion  reso- 
nances, and  for  the  ion-ion  hybrid  resonance  are  presented. 

PACS  Index:  52.25.Ps,  52.35.Hr,  52.50.G 

Permanent  address:   College  of  Staten  Island,  City  University  of  New  York, 
St. George  Campus,  Staten  Island,  NY  10301. 
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I.  INTRODUCTION 

Heating  of  plasmas  by  waves  at  freqaencies  matchiag  some  ion  cyclotron 
harmonic  resonance  has  attracted  attention  both  experimentally  and  theoreti- 
cally. In  particular,  the  second  harmonic  resonance  heating,  with  or  with- 
out minority  ion  species  at  the  fundamental  resonance,  has  been  used  experi- 
mentally with  success  to  heat  plasmas.  Using  frequencies  matching  a  higher 
ion  cyclotron  harmonic  resonance  may  be  attractive  since  smaller  antennas 
for  launching  the  waves  are  appropriate.  Successful  use  of  these  heating 
schemes  depends  on  the  energy  absorption  efficiency,  as  well  as  the  accessi- 
bility of  the  incident  waves  to  the  resonance  region.  The  purpose  of  this 
paper  is  to  present  a  study  of  propagation,  absorption,  and  mode  conversion 
of  waves  across  cyclotron  harmonic  resonances  for  a  perpendicularly  strati- 
fied medium  for  waves  propagating  nearly  perpendicularly.  When  the  waves 
propagate  nearly  perpendicular  to  the  equilibrium  magnetic  field,  the  stan- 
dard geometrical  optics  approximation  is  expected  to  breakdown  near  the  res- 
onance layer.  Associated  with  this  failure,  reflection  and  mode  coupling  ef- 
fects may  also  become  significant,  so  that  a  full  wave  analysis  is  required. 
In  Refs.1-4  we  have  studied  similar  problems  for  the  electron  cyclotron  har- 
monics. By  use  of  a  boundary  layer  expansion  about  the  resonance  surface  we 
were  able  to  obtain  full  wave  solutions. 

The  work  here  overlaps  some  of  the  extensive  analysis  of  Swanson,*'*  and 
Colestock  and  Kashuba',  but  it  differs  in  several  aspects.  We  prefer  to 
treat  the  system  by  boundary  layer  expansions  rather  than  inference  of  wave 
equations  from  uniform  medium  dispersion  relations  and  energy  conservation 
relations.  We  believe  that  the  boundary  layer  expansion  technique  gives  the 
cleanest  and  simplest  set  of  approximate  equations  for  the  study  of  wave 
propagation  near  a  cyclotron  resonance.  By  evaluating  the  current  from  the 
Vlasov   equation  with  our  method,   we  obtain  directly  the  wave  equation   for 
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the  electromagnetic  wave  field,  and  not  more  abstract  wave  equations.  Since 
we  are  finally  interested  in  complete  calculations  of  full  wave  fields,  see 
e.g.  Refs  1-3,  rather  than  less  complete  but  far  simpler  variational  calcu- 
lations, we  focus  on  a  full  wave  treatment.  By  the  use  of  full  wave  equa- 
tions, we  obtain  a  group  of  exact  properties  of  the  solutions.  Specifical- 
ly, we  obtain  properties  that  are  independent  of  the  equilibrium  distribu- 
tion function,  provided  the  distribution  function  is  isotropic  in  velocity 
space.  We  obtain  reciprocity  theorems  and  we  find  general  connections  be- 
tween the  full  wave  solution  to  a  scattering  problem,  and  the  geometrical 
optics  solution  to  the  same  problem. 

Recently,  one  of  us  (K.I.)  has  investigated  the  wave  propagation  across 
the  higher  (i.e.,  third  and  higher)  electron  cyclotron  harmonic  layers  using 
the  same  method.*  It  was  shown  that  the  relativistic  corrections  wash  out 
resonance  effects,  and  geometrical  optics  approach  becomes  applicable.  In 
that  case,  to  the  lowest  significant  order  in  the  expansion  parameter,  the 
ratio  of  the  electron  thermal  speed  to  the  speed  of  light,  we  found  that  the 
cold  plasma  approximation  holds,  and  the  effects  of  the  higher  harmonic  res- 
onance are  negligible.  In  the  present  work  we  address  ourselves  to  similar 
questions  in  the  case  of  ion  resonances.  By  making  use  of  a  similar  boundary 
layer  analysis,  we  show  that,  unlike  the  electron  case,  there  exists  certain 
parameter  domains  where  the  geometrical  optics  approximation  fails  and  fall 
wave  treatments  are  needed.  For  the  case  of  ions  we  expect  relativistic 
corrections  to  be  negligible,  but  as  a  matter  of  completeness  we  determine 
the  extremely  small  range  in  nn  in  which  they  might  appear. 

We  consider  in  detail  three  essentially  distinct,  but  closely  related 
problems.  Nth  harmonic  resonance  N>_2,  fundamental  cyclotron  resonance,  and 
ion-ion  hybrid  resonance  with  a  majority  ion  species  at  second  harmonic  res- 
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onance  and  minority  species  at  fundamental  resonance.  We  obtain  wave  equa- 
tions appropriate  to  each  case.  In  all  these  problems  we  have  a  fast  mode 
—  an  extraordinary  mode  for  N>^2  and  a  magnetosonic  mode  for  N=l  —  which 
propagates  on  both  sides  of  resonance,  and  a  slow  mode  —  a  Bernstein  mode 
for  N>_2  —  which  propagates  on  only  one  side  of  resonance.*  We  characterize 
the  wave  propagation  problem  in  terms  of  coefficients  of  wave  transmission, 
mode  converted  transmission,  reflection,  and  mode  converted  reflection.  We 
prove  a  series  of  exact  results  corresponding  to  the  generalization  of  the 
reciprocity  theorems  of  Chin  and  Man.*  In  another  vein,  for  the  cases  of 
second  harmonic  resonance,  fundamental  resonance,  and  ion-ion  hybrid  reso- 
nance we  obtain  exact  formulas  for  the  transmission  coefficient  for  incident 
fast  modes,  and  show  that  for  a  fast  mode  incident  from  the  high  field  side 
no  fast  mode  is  reflected.  We  believe  that  our  method  of  proof  of  these 
latter  results  would  also  extend  directly  to  higher  harmonic  resonances. 
Additionally,  for  purely  perpendicular  propagation  for  all  cases  except  ion- 
ion  hybrid  resonance  we  obtain  integral  representations  of  the  solutions. 
Finally,  we  examine  these  cases  numerically  to  obtain  information  on  plasmas 
of  some  laboratory  interest. 

The  organization  of  this  paper  is  as  follows.  Section  II  contains  the 
formulation  of  the  boundary  layer  analysis,  and  the  derivation  of  the  cur- 
rent density  from  the  Vlasov  equation  for  Nth  ion  cyclotron  harmonic  reso- 
nance. We  use  the  current  density  with  Maxwell's  equations  to  derive  field 
equations  valid  within  the  boundary  layer  in  Sec. Ill  for  the  case  of  Nth 
harmonic  resonance  N>,2.  We  also  establish  the  connection  of  the  solutions  of 
these  equations  with  those  given  by  geometrical  optics.  Furthermore,  the 
wave  equations  enables  us  to  prove  the  reciprocity  relations  in  general 
terms.  These  relations  state  that  (i)  the  transmission  coefficient  for  the 
X-mode  is  independent  of  the  direction  of  incidence,  (ii)  the  reflected  mode 
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coupling  coefficients  are  the  same  in  the  cases  of  fast  mode  and  the  slow 
mode  incidences.  and  (iii)  the  transmitted  mode  conversion  coefficients  are 
the  same  in  the  cases  of  fast  and  slow  mode  incidences.  The  validity  of 
these  relations  does  not  depend  explicitly  either  on  parallel  index  of  re- 
fraction or  on  the  equilibrium  distribution  function.  The  method  of  proof  of 
Chiu  and  Mau*  which  leads  to  similar  reciprocity  relations  cannot  be  direct- 
ly applied  in  our  case  as  it  exploits  the  symmetry  of  the  full  system  which 
may  not  be  retained  in  the  approximate  equations  derived  for  the  resonance 
layer.  Two  other  results  found  in  the  special  case  of  zero  parallel  refrac- 
tive index  are  that  the  reflection  coefficient  is  zero  for  a  high  field  side 
incident  fast  mode,  and  that  the  transmission  coefficient  for  this  mode  is 
the  same  as  that  predicted  by  geometrical  optics.  We  also  show  that  these 
results  hold  generally  for  second  harmonic  resonance.  We  suspect  the  re- 
sults are  true  for  all  harmonics,  and  the  proof  would  be  an  extension  of  the 
theorem  in  the  Appendix. 

In  Sec. IV,  we  study  the  problem  of  fundamental  resonance  with  one  ion 
species.  We  show  that  there  a  thin  boundary  layer  about  the  fundamental 
resonance  layer  in  which  geometrical  optics  fails.  The  dispersion  relation 
in  a  uniform  medium  predicts  a  bifurcation  point  within  this  resonance  lay- 
er, where  the  fast  and  slow  branches  of  the  magnetosonic  mode  coalesce. 
Partially  on  the  basis  of  the  cold  plasma  model,  one  expects  that  a  fast 
mode  is  largely  unaffected  by  the  fundamental  cyclotron  resonance,  and  that 
the  slow  mode  may  be  absorbed.  The  Vlasov  model  is  more  complex,  with  mode 
conversion  and  tunnelling  occurring.  Nonetheless,  this  crude  picture  per- 
sists. The  mode  coupling  rapidly  vanishes  outside  a  narrow  cone  of  incident 
waves.  The  fast  mode  incident  from  the  low  field  side  is  also  partially 
reflected.    The   proof  of  the  reciprocity  relations  for  this  case   is   also 
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glven.    The   fast  mode  incident  transmission  coefficient  is  given   exactly, 
and  the  high  field  side  incident  fast  mode  has  no  reflected  fast  mode. 

It  is  well-known  that  the  absorption  of  waves  propagating  across  the 
second  harmonic  layer  is  significantly  improved  by  the  inclusion  of  a  small 
amonnt  of  another  ion  species  at  the  fundamental  resonance  within  the  same 
resonance  surface.  In  Sec.V,  we  study  this  problem  by  our  boundary  layer 
technique.  Clearly,  when  the  transmission  coefficient  is  near  unity,  the 
absorption  is  necessarily  poor.  We  show  that  by  inclusion  of  minority  ions 
we  can  significantly  reduce  the  transmission  coefficient.  A  substantial 
fraction  of  the  incoming  wave  energy  is  then,  in  fact,  absorbed  by  ions. 
The  addition  of  the  minority  species  even  farther  reduces  the  parameter 
range  in  which  geometrical  optics  can  be  applied  across  resonance.  In  fact 
in  some  cases  the  reflected  mode  coupling  remains  significantly  large  even 
when  the  parallel  refractive  index  is  Increased  close  to  its  left-hand  cut- 
off value.  It  is  largely  because  of  this  property  that  the  absorption  of  the 
waves  Incident  from  the  low  field  side  is  higher  than  that  of  the  high  field 
side.  Again,  reciprocity  theorems  and  exact  values  of  transmission  coeffi- 
cient are  given.  In  Sec. VI,  we  present  our  numerical  results  for  the  var- 
ious cases  studied  in  this  paper.  We  show  that  the  third  harmonic  resonance 
can  be  used  to  heat  tokamak-like  plasmas  by  low  field  side  Incident  fast 
waves.  The  high  field  side  Incident  fast  (X)  or  slow  (Bernstein)  modes  may 
provide  acceptable  absorption  rates  if  the  parallel  refractive  index  is  suf- 
ficiently large.  Absorption  of  fast  magnetosonic  waves  at  the  fundamental 
resonance  could  be  significant  for  dense  and/or  large  plasmas,  whereas  the 
slow  mode  is  readily  absorbed  even  when  n,,  is  quite  small.  We  study  the 
parameter  range  corresponding  to  PLT,  JET,  and  CIT  devices  for  the  case  of 
ion-ion  hybrid  resonance.  We  find  that  when  the  minority  number  density  is 
approximately  1%  of  the  electron  number  density  the  low  field  side   incident 
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fast  mode  absorption  rate  peaks  strongly.  When  nii  is  sufficiently  small 
most  of  the  absorbed  incident  energy  is  deposited  to  the  minority  species. 
We  find  that  in  the  case  of  CIT,  absorption  of  the  high  field  side  inciden* 
fast  or  slow  modes  is  qaite  negligible. 

Section  VII  has  a  brief  discussion  of  the  results  and  concluding  re- 
marks. An  appendix  contains  a  useful  theorem  on  as3niiptotic  expansions  of  a 
solution  of  ordinary  differential  equations. 

II.  EVALDATION  OF  CURRENT  DENSITY 

A.  The  Vlasov  Equation  within  the  boundary  layer 

Within  the  ion  cyclotron  resonance  layers  we  represent  the  plasma  ions 
by  the  linearized  Vlasov  equation,  and  we  solve  the  Vlasov  equation  by  per- 
turbation expansion.  The  consistency  of  the  expansion  is  not  fully  demon- 
strated until  the  next  section.  The  equilibrium  magnetic  field  B  is  as- 
sumed to  be  along  the  z  direction,  and  its  magnitude  varies  slowly  along  the 
X  direction.  The  time  and  z  dependences  are  assumed  to  be  harmonic  in  the 
form  exp[-i(i)(t  -  znnco/c)].  We  normalize  the  velocity  variable  with  respect 
to  the  ion  thermal  speed  v  ,  and  introduce 


■  =  v.  w, 

3  3 

gdw=evf.du,  (1) 

3        3 
n.  Gdw=F.  du, 
1  1 

where  f.  and  F.  denote  perturbed,  and  equilibrium  ion  distribution  func- 
tions,  respectively.  We  assume  that  G  is  isotropic  in  the  velocity  space. 
Within  the  Nth  harmonic  resonance  layer  we  introduce  the  stretched  spatial 
variable 

1  1  _  QU/LK 
^   8  ^N      u)  ^'  ^^' 
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where   the  value  of  5  will  be  determined  to  assure  consisteocy  of  the  expan- 
sion, and  u)  =  NQ(0).   We  write  the  ion  Vlasov  equation  in  the  form 


(1  -i.8  Aj)gj  .X(gj^^  .gj_^)  =b.. 


(3) 


where   we  have  expanded  g  in  a  Fourier  series  in  gyrophase  angle  0,   and   we 
set 

w  =  (W|  cos  i,    wi  sin  i.    Wn  ), 

A  jt  ff~l  J-      2..-1   w 

A.    =   J5   -   1,8        n||W,|     +  n   8      —   . 

b.    =  ^  YX.[i-  ^    (5:4      +   8^4,)    +   f^     8°   E,|  ]. 
j         4ni    '    1^2    3w|      J      -  j      +  aw.,  "j       ||  " 

Q   =  Z.eBQ(x/L)/m^c, 
r  =   (1   +  .^  w^)^^^ 

n  =  Vj/c, 

E     =  E     T      IE   ,    Em     =  E    , 
T  X  y        II  z 

where  n^  is  the  cold  plasma  refractive  index  evaluated  at  the  resonance.   We 

specify  n   in  the  next  section.    In  order  to  evaluate  the  current  in   terms 
c 

of  the  field  variables,  we  must  calculate  g_,>  J  =  0,7l  from  the  system  (3). 
In  the  present  work  we  are  only  interested  in  the  lowest  significant  order 
contributions  of  the  resonance  effects.  Thus  we  keep  only  the  terms  of 
order  zero,  namely 

gj°^  =  gj  {1  +  0(e^)}.  (5) 

and  suppress  the  superscript  (0)  for  brevity,  so  that  (12)  can  be  written  as 


1 
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with  K(55)  replaced  by  K(0)  in  L. 

The  right  hand  side  of  the  equation  corresponding  to  the  resonant  com- 
ponent g^  is  zero  except  for  the  case  N=l.  Thus  we  must  treat  the  case  of 
the  fundamental  resonance  separately.  We  first  consider  the  case  when  N>^2. 
B.  Ion  Current  when  N  >_  2 

We  avoid  making  explicit  assumptions  about  the  relative  orders  of  the 
Fourier  components  (the  ordering  will  become  evident  at  the  end  of  our  anal- 
ysis), except  that  we  hjrpothesize  that  g  tends  to  zero  for  large  |j|.  By 
induction  we  obtain  the  relation 

8_2  =  -i^^8-l  ^0(g_„).  M>3.  (7) 

from  which  it  follows  that  we  may  drop  the  term  involving  g  in  the  equa- 
tion corresponding  to  j  =  -1  in  (6).  Next,  we  express  g  in  terms  of  g  and 
gj.,  by  making  use  of  (6)  for  j=2,3  ,  •  •  •  ,N-1 ,  yielding 

N  r   ^  (-1)'^  mN-2  r^-1  ,«, 

We  may  relate  g„  to  g  by 

We  complete  our  truncation  process  by  expressing  g-vr_.  in  terms  of  g^^.  We 
accomplish  this  by  using  the  relation 

»m^N  =  I  K.K-1   ^  0(^Vn.M>'  ">1'  »'>°'  (^°> 

which  again  may  be  established  by  induction.   For  m  =   N-1 ,  and  for  M  snffi- 

M  7 

ciently  large  so  that  e  g2N+i+M  ~  0(c),  the  latter  gives 
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N^-^   ,N-1 

We  now  obtain  the  solation  of  the  Vlasov  equation  by  asing  the  relations 
(7),  (8),  (9),  and  (11),  By  making  use  of  (6)  with  j=0,  j=-l,  and  j=l,  we 
express  g   in  terms  of  g  ,  and  thus  in  terms  of  g^^,  giving 

N-1 
(-N)   ^.N-1    ^  -  ,_, 

where 

We  next  substitute  (12)  into  (9),  and  use  (11)  to  obtain 

N-2 
-  .  (-N)*"  ^    .N-1  - 

^  ^N  8n  =  -  (N-2)!   ^    ^-  ^^*^ 

We  have  thus  expressed  the  resonant  component  in  terms  of  the  field  varia- 
bles alone,  and  thus  we  can  evaluate  every  g  in  terms  of  b,.  In  particu- 
lar, 

C  ..      N^^"^      /-N-1     1     .N-1  -  ,,., 

*1  =  ''  "■  (N-l)l  (N-2)  I   ^      ^   ^H 

The  second  term  in  the  right-hand  side  of  (IS)  represents  the  contribution 
of  the  Nth  harmonic  resonance  to  the  lowest  significant  order  in  the  expan- 
sion parameter.  Similar  relations  for  g  and  g_.  are  also  obtained.  The 
corrections  are  of  higher  order  in  e    in  these  components: 

N   r  .      N   /•  ,        N^^^    .N   1    ,N-1 


h   =  \-  ii^^^l-  N^^^-l   ^   j(N-l 


^N-2 

)I1''     8  A^       ^         (1(f) 


2  ~        ~    ^' 


N  /v      f  u   ^   N   fZ  V  X  ^      N^^      /-N+l    1    M-1  . 
8-1  =  Wn^^-1   -  -^  ''O  ^  N=r  ^  ^^  ""  (N+D!  (N-1)  I  ^     R  A        1* 

We  obtain  for  the  ion  current  density  from 
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T-ion  rd   w  _ion    rd  w  ,,„^ 

J^  =  /—  wj^g^j.   J||  =   /V'll'o'  (1^) 


so  that 
4ni 


0.  ^11  =  "^i^l 


.  2ic^  X.  B  D^a^_^^/-lE^  -  ie^^^/^^E,  -  iea^.^^^D^E,,  ). 


inl  T     JL  T  ^i?   JL  1  2,  2.n2c  X 


^12  _N+1,     _N-1_   ^12  N(N-t-2)       „N+1_,     .      „N„   , 
-  Ck,  X,  -c   D    (a„  „D   E^  +  -e   -^^ fa^,^,  „D    E  -  iea„  n  E„  ). 


4ni  -      _N_  T   /c      N  1  2  .  2^_2„  , 

where  <•••>  denotes  velocity  space  average  with  respect  to  the  eqailibriom 
distribntion  function,  and 


(18) 


5 

(19) 


"^  [(N-l)!]2 


D  =  K(0)77  . 

The  functions  a    are  defined  by 
n.m 

1^  f  d  w   2n+l   m  30  ,,„, 

"n.m  =  -  ?^  J  "i;;  'i     "11  j;;^-  (^o) 

We  now  retarn  to  the  case  omitted  in  this  section,  N  =  1. 
C.  Fundamental  ion  resonanco 

The  fundamental  ion  cyclotron  resonance  problem  differs  significantly 
from  the  other  harmonic  resonances  which  we  have  just  studied.  The  basic 
difference  stems  from  the  fact  that  now  the  resonant  cold  plasma  dielectric 
tensor  element  e  becomes  unbounded  as  one  approaches  the  fundamental  reso- 
nance.  Thus,  it  becomes  possible  that  E  decreases  near  resonance,  and  cor- 
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respondinglj  the  sole  effect  of  resonance  may  be  to  shield  oat  E  .  Addi- 
tionally, although  no  ion  Bernstein  modes  propagate  near  the  fundamental 
resonance  layer,  other  modes  do  exist,  and  it  is  necessary  to  investigate 
carefully  the  uniform  medium  dispersion  relation  in  detail  near  resonance  in 
order  to  discover  what  other  modes  are  present.  We  give  this  analysis  as 
part  of  Sec. IV. 

We   return  to  the  system  (l)-(4),  with  N  =  1,  and  we  take  n   to  be   the 

c 

cold  plasma  refractive  index  evaluated  at  the  resonance,  which  we  specify 
later.  Again,  we  solve  the  chain  equations  (3)  order  by  order  in  terms  of 
field  variables.   To  the  lowest  significant  order  in  e  we  obtain 


«-2 


3  -^^-l' 


gj  =  ^Si' 


so  that 

8Ag   =b   -Xb^  2  ^Vi- 

\  (21) 

'O  =  ''O  "  2  -^  ^-1  "  -^  'l' 

lu      l/-u^lf2 

For  any  cyclotron  resonance  but  fundamental  we  could  determine  the  relation 
between  5  and  e  by  demanding  that  current  densities  of  order  one;  now  we 
must  turn  to  Maxwell's  equations  to  find  the  relation.  We  calculate  the  ion 
current  density  without  making  any  assumption  on  parameters  and  on  the  rela- 
tive ordering  of  the  variables,  other  than  requiring  that  (21)  be  correct  to 
leading  order.  The  components  of  the  ion  current  density  are  then  obtained 
using  (18) .   We  find 

2 

^^  J^  =  -7   X^n  «E  +  fr  X.a,  ^D^E  +  if  X.a,  ,DE„  , 
u)    +    5   1  0,0  +    45   1  1,0    -    8   1  1,1   II  .   . 

2  4  2 

~  ^-  =  °  «  ^i"i.o^  ^  ^<-r  -  ^ ""  '^2.0°  ^^-  -  JT  ^°  'i.i^^ii  • 
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2  2 

*^  J"  II     =  T7  XDa«   ,E     -   if-  X.Da,    ,D^E     +  X.(-l    +   fr  Da„   ,D)E,|   , 
CO        II         26      1     0.1   +        85        1      1,1        -  1  25        0,2         II 

where  a     is  as  given  in  (20)  with  N  =  1. 

n,m  ° 

III.    HARMONIC    RESONANCE    WITH    N    >    2 

A.    Field    Eqaationt    within    the    Resonance    Layer 

We  calculate  electron  corrent  evaluated  at  the  resonance  from  the 
cold  plasma  conductivity  and  the  additional  assumption  of  charge  neutrality 
of   the    system,    i.e. 

^J   '^*^*'   =  ±NX.  (23) 

(0      7  1 

and   substitute   it  with  the  ion  current  density  (18)  into   Maxwell's   equa- 
tions.  We  normalize  the  problem  by  assuming  that  E  =0(1),  and  find  that  the 

2N-2 
stretching   parameter  5  must  be  chosen  to  be  e     so   that   the   resonance 

effects  will  be  present  at  the  lowest  significant  level  in  the  boundary  lay- 
er expansion.   When  we  drop  the  terms  of  order  m  /m. ,  and  eliminate  E   with- 

e   1  — 

in   the  resulting  equations  we  obtain  the  following  equation  for  E  which  is 
valid  in  the  Nth  harmonic  ion  resonance  layer 

2         2      N-1     N-1 
[1  +  d''  +  (1  +  d''  +  n)d   ^  A  d   ^]E^  =  0,  (24) 

where  d  =  d/ds,  and 

n  (0 

^  =  5  f(f)  77  =  -5%s^/'=' 

pi 

\  =  (E__  -  4 )(%+  -  4  ^^^Sx  -  4  ^'  Sx  =  2^'++  ^  '-^' 


^  .  .  N^ 


2  2 

]i  =   (e__  -  nji  )/(e^^  -  n||  )  , 

X.               ^   n^^N-l    ,    ,,N+1,.,,-,2N-1 
.   ^  1  r c]  (-1)        (N/2) 

I  2    ,    L  X.J  ,,-,  ,,,,2  N-1.0    • 

(e        "  °||  )  ^  I  (N-1)  !] 
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The  components  e,,  given  above  correspond  to  the  cold  plasma  dielectric  ten- 
sor elements  evaluated  at  the  resonance.   and  n   is  the  Appleton-Hartree  in- 

c 

2 
dex  of  refraction  again  evaluated  at  the  resonance.   We  note  that  if  nii  << 

X,  then  we  can  drop  nii  in  the  definitions  given  in  (25).   However,   we  see 

from  (4)  that  even  when  these  corrections  may  be  dropped,  when 

n||  =  0(8/t,), 

the   effects  of  nn  on  wave  propagation  are  significant.    It  is   remarkable 

that  the  relativistic  corrections  are  also  significant  within  the  layer  pro- 

2 
vided   that  the  more  stringent  condition  hold  that  t)  /5  be  of  order  anity  or 

less.  However,  these  corrections  are  not  large  enough  to  suppress  the  reso- 
nance effects  as  they  do  in  the  case  of  higher  electron  harmonics.  as  was 
demonstrated  in  Ref.4.  Hence  the  relativistic  effects  are  significant  with- 
in a  thin  propagation  cone  if  the  following  ordering  holds: 

X.  =  0(n   ^~^).      and   nji  =  0(ii)  . 

To  show  that  the  asymptotic  solutions  of  Eq.(24)  can  be  connected  to 
the  geometrical  optics  solutions,  we  next  derive  the  asymptotic  expansions 
of  the  solutions. 

B.  Asymptotic  Expansion  of  the  Solutions  of  Inner  Equation 
For  1^1  large,  we  can  write 

1   .   2N-2.  ,,,. 

where  <•••>  denotes  velocity  space  average  with  respect  to  the  equilibrium 
distribution  function  G.  so  that  Eq.(24)  becomes 

{1  +  d^  -  [b  +  a(l  +  d^)]d^"^  -  d^"^}E^  =  0.  (27) 

where 
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b  =  (la, 

X.       _  n^,N-l  ,  ,vN+l,.,,-,2N-l     ,„  ,  (28) 

a  =  ^-^     -^     ^-'>        ^^^'\  <.'^-h    (n^  5  L  J,  . 

(e   -  nt,  )  ^   ^1  [(N-l)!]^      -L        '^  *" 

zz     II 

We  note  that  (27)  also  represents  the  case  of  purely  perpendicolar  propaga- 
tion in  the  nonrelativistic  limit.  For  X.  large,  the  case  of  interest,  |i<0, 
so  that  a  and  b  are  of  opposite  sign,  and 

a  =  lal  (-1)^*^.  and  b  =  |b|  (-1)^.  (29) 

We  can  write  the  solution  of  (27)  in  the  form 

u  =  -/  e^'   p^-1  (a  +  -^)  e^^P^  dp,  (30) 

1+P 

where  the  contour  of  integration  in  (30)  is  to  be  specified  later,  and 
2N-1       N-2     .   2j+l     ,,  , 

P(p,  =  -  .  ^   -  IM^^^(-uJ  ^    .  M  ,,,^ 

The  integrand  of  (30)  tends  to  zero  as  |p|  — >  <»  within  the  sectors 

2k-l  n  ,       ,  2k+l  n   ,    .  ^„   ..  .  . 
W^   2  ^  "«  P  ^  2N:T  2'  ^"  ^^^   °"  integer. 

The  sectors  corresponding  to  k  =  7(N-1)  are  of  particular   importance.    The 

contours   extending  to  infinity  within  these  sectors  will  be  seen  to  include 

the   contributions  of  the  propagating  ion  Bernstein  waves.    For  1st   large, 

only   the  contributions  of  the  branch  points  at  p   =  ri,  and  the   stationary 

b 

points  p  need  be  considered.   The  stationary  points  are 


P.  =  (f)^^  ^  z.  ,   k=0,l,---,2N-l,  (31) 

s    a       k 


where 


z.    =  exp(in  JTT)  .    for   s/a   >>    1, 

2k-l 
=  exp(in  TjZr) .    for   s/a   <<   -1. 
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The  contribations  of  the  stationary  points  generate  either  asymptotically 
small  or  asymptotically  large  solutions  except  for  the  branches  p  =  7ilp  I, 
which  occurs  only  for  s>0  (high  field  side).  We  find  that  there  are  only 
three  (out  of  2N)  linearly  independent  solutions  which  are  physically  rele- 
vant for  the  scattering  problems.  The  contours  C  ,  C  ,  and  C,  which  gener- 
ate these  solutions  are  shown  in  Fig.l.  The  other  linearly  independent  so- 
lutions all  are  asymptotically  unbounded  as  s  — >  +">  or  s  — >  -".  Moreover, 
these  solutions  are  linearly  independent. 

By  using  standard  methods  we  find  that  the  asymptotic  expansion  of  the 
three  linearly  independent  solutions  for  |s|  large  are 

lb|5n/4  ^      ™„ 
u  =  p   e        +  <l_  t  HFS, 

lb|3n/4   ,„  ^"' 

=  p   e       ,   LFS, 


and 


|bl3n/4.  Ibln/2   „-lblH/2     +  „  ^  +  „   +  „   HFS 
u  =  e       (e       -  0       )  (p   +  p_)  +  q  +  q_,  tn<», 

^  (33) 

=  0.  LFS, 


where   HFS  and  LFS  stand  for  high-field  side  (s>0)  and  low-field  side  (s<0), 
respectively,  and 


p^  =  T  I  n±ilbl/2)    exp[Ti(nj  +   s^   +  ^logZls^l)!. 

'l   =   '^  -  iS^l  -    ''''    ji^jTI' 

s^   =   s  -    lal    -   (N-l)lbl. 

1/2 
Iff    a  I  lb  Iff  r     -f     I        i2M-2  Tf  /tKt^t  \^^ 


(34) 


We  have 


,2    2|b|ff,  lblff/2    -lblff/2,   *Ps' 


Iq^/p^j   =  e     (e       -  e       )  ^^^_^^    .  (35) 
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In   closing  this  section  we  note  that  since  q    represent  incoming  and   out- 
going  ion  Bernstein  modes  we  identify  the  eikonal  function  for  this  mode  as 

2N— 2 
sip  I  „_  .    By   calculating  the  gradient  of  this  function  we  find  that   the 

refractive  index  n^  of  this  mode  is  related  to  n  by 
Q  c   ' 


n./n   =  Ip  I  =  |s/a|^^  ^  .  (36) 

Q  c      s 

One  can  easily  re-derive  (36)  by  making  use  of  uniform  medium  dispersion  re- 
lation evaluated  near  the  Nth  resonance  layer. 

In  order  to  evaluate  the  mode  coupling  coefficients  we  need  to  know  the 
ratio  of  the  components  of  the  energy  flux  vectors  along  the  direction  of 
propagation  for  the  fast  and  slow  modes.  The  energy  flux  for  the  fast  mode 
can  readily  be  obtained  from  the  Poynting  theorem,  giving 

F^  =  7-  Re(B*  X  B)   =  ^  n  (e   /e   )^|E^|^.  (37) 

X   4n  X   4n   c   XX   —    + 

The   energy  flux  associated  with  the  slow  (Bernstein)  mode  is  derived  as   in 
Ref.l  with  ions  replacing  electrons,  and  we  find 

Fo  =  r—  «    lEj^  (38) 

Q   4n   n   xx   + 

It  is  clear  from  the  above  relation  that  for  the  ion  Bernstein  mode  with   X. 

1 

>>1   there  is  a  180  degree  phase  shift  between  the  direction  of  propagation 
and  the  direction  of  energy  flow,   since  e    is  negative  near  the  resonance. 

XX 

We  readily  find  that 

,Q|2 


^X     'Ps'  Ie^I 


We  now  calculate  the  solutions  of  basic  scattering  problems  for  this  special 
case  of  the  perpendicular  propagation  by  taking  appropriate  linear  combina- 
tions of  the  three  solutions  (32),  (33).   We  readily  find  that  when  the  fast 
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mode  is  incident  from  the  high-field  side  the  transmission  coefficient  for 
the  energy  flax  is 

T=e-"'^'.  (40) 

We  also  find  that  the  mode  converted  energy  flax  reflection  coefficient  R^  = 
1  -  T,  so  that  there  is  no  reflected  X  mode  in  this  case.  The  absorption 
within  the  resonance  layer  is  zero.  If  the  fast  mode  is  incident  from  the 
low-field  side  we  find  that  the  transmission  coefficient  for  the  energy  flax 

is  again  given  by  (40).    There  is  a  reflected  fast  wave  in  this  case  and  IL 

2 
=  (1  -  T)  .    The  transmitted  mode  converted  coefficient  for  the  energy  flax 

T  =  T(l-T),  and  the  absorption  is  zero. 

When  a  Bernstein  mode  is  incident,  a  Bernstein  mode  is  reflected  with 
reflected  energy  flax  coefficient  R^  =  i.  where  T  is  as  defined  in  (40).  We 
also  find  that  the  mode  converted  transmission  coefficient  and  the  mode  con- 
verted reflection  coefficient  are  the  same  as  those  jast  given.  Again  there 
is  no  absorption  within  the  resonance  layer. 

As  we  have  jast  shown,  there  is  no  net  absorption  of  the  incident  wave 
energy  by  the  plasma  within  the  resonance  layer  when  the  relativistic  cor- 
rections are  dropped  and  the  propagation  is  parely  perpendicalar.  Althoagh 
the  energy  in  Bernstein  waves  generated  by  mode  conversion  may  be  eventual- 
ly absorbed  by  the  plasma,  the  energy  deposition  dae  to  such  processes  does 
not  occar  in  the  resonance  layer,  and  hence  we  do  not  include  it  here  as 
energy  flax  absorption.  Without  the  extremely  small  relativistic  correc- 
tions we  obtain  energy  deposition  within  the  resonance  layer  only  with  the 
inclusion  of  a  small  parallel  refractive  index.  Unfortunately,  unlike  the 
special  case  just  studied,  it  is  not  possible  to  construct  the  solution  gen- 
erally in  terms  of  quadratures.    It  is  therefore  necessary  to  perform  a  nu- 
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merical  study  for  this  case  similar  to  that  of  Refs.   1-2.    We  present  this 
study  for  N=3  in  Sec. VI. 

In  order  to  examine  (24)  more  efficiently,   we  write  it  as  the  system 

T'  =  M  •  T.  (41) 

where  T  =  (y-i 'y^ » *  *  * 'y^N^ '   *°^  ^1  ~  °'   ^®  matrix  M  has  only  the  following 
nonzero  elements: 


M.     =  1,  for  j=l.---,N-2.  or  j=N, • • • .2N-3 ,  or  j=2N-l. 


V,^^=l/A(s) 


**2N-2.1  =  "*'2N-2.2N-1  =  ^'  ^^2) 

Vl  =  ^' 
*'2N.2N-1  =  -^l*^^'- 

It  was  shown  in  Ref.  1  that  a  Riccati  transformation  is  useful  not  only  to 
reduce  the  order  of  the  differential  equation,  but  also  to  study  the  asymp- 
totic behavior  of  the  components  of  T.   We  define 


Zj  =  yj+i/yj.  j=i.---.N-i.  (43) 

and  derive  the  nonlinear  reduced  system  of  differential  equations  for  z,. 

T.\    =  z+i  -  Zj  z  ,  j=l,---.N-l,  and  j=N,  •  •  •  ,2N-2,  and  j=2N-2, 

_  i 

*N-2   A  ^N  "  ^1  ^N-2'  ,,,, 

144; 

*2N-3  "  ^  "  ^2N-2  "  ^1  ^2N-3' 
Sn-1  =  ^^  -  <l'-»^^^2N-2  "  ^1  ^2N-1- 

Since  z'  should  tend  to  zero  for  |s|  large,  we  readily  obtain  for  |s|  large 

r.  =  z^,  j=2,«-',N-2, 

z  =  A  zj,  j=N-l,---,2N-3, 
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i  i   2N-2  (.,. 

1  -  ^2N-2  =  ^  "l    '  ^*^^ 

^2N-1  "  ^1  ^2N-2' 

We   see  that  if  z      is  of  order  one,   then  z,^,  =  !•   which  in  turn   implies 
that  z   =  -1,   and  z,»t_,  =  z,  .   Clearly  we  have  obtained  the  X  mode.   Alter- 

natively,  if  z     =  1/A,  then  Z2N-2  ^  ^^^i'    ^^^   ^2N-1  ^  ^^^i'      ^^  selecting 

-1/ (2N-2) 
any   one   of  the  2N-2  branches  of  A  we  can   obtain   the   asymptotic 

values  of  the  remaining  solutions 


C.  Reciprocity  Relations 

The  reciprocity  relations  which  we  have  obtained  in  the  special  case  of 
purely  perpendicular  propagation  in  the  nonrelativistic  limit  are  also  true 
generally.  These  relations  state  that  (i)  the  transmission  coefficient  is 
the  same  when  the  fast  mode  is  incident  from  the  low  or  high  field  side, 
(ii)  the  reflected  mode  conversion  coefficients  are  the  same  in  the  case  of 
fast  and  slow  mode  incidences,  and  (iii)  the  transmitted  mode  conversion 
coefficients  are  also  the  same  in  the  fast  or  slow  mode  incidences.  In  addi- 
tion to  these  reciprocity  relations  there  are  two  other  properties  which  we 
have  proved  in  a  special  case,  which  state  that  for  X  mode  incident  from  the 
high  field  side  there  is  no  X  mode  reflected  and  the  X  mode  transmission 
coefficient  is  exactly  the  geometrical  optics  result  (41).  In  the  Appendix 
we  present  a  theorem  which  includes  these  results  as  a  special  case  for  N=2, 
for  essentially  arbitrary  parallel  refractive  index  and  distribution  func- 
tion. The  method  of  the  Appendix  could  almost  assuredly  be  applied  for  ar- 
bitrary N>2  to  obtain  the  same  conclusion. 

Both  the  reciprocity  theorems  and  the  exact  values  of  the  appropriate 
reflection  and  transmission  coefficients  were  found  in  the  numerical  calcu- 
lations of  second  harmonic  electron  resonance  problem  in  Ref.2.   With  slight 
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modif ications,   the   theorem  of  the  Appendix  and  the  reciprocity  proofs   of 
this  section  coald  also  be  applied  to  electrons. 

In  Ref.3>  for  the  case  of  electron  fundamental  resonance,  we  showed  the 
equality  of  the  two  transmission  coefficients  by  using  the  Wronskian  rela- 
tion. Since  the  order  of  the  differential  equation  in  that  case  is  two,  the 
Wronskian  relation  provides  a  bilinear  relation  between  any  two  solutions  of 
the  system.  In  the  case  of  the  Nth  resonance,  the  Wronskian  relation  is  not 
very  useful,  since  it  relates  all  2N  linearly  independent  solutions  of  the 
system.  In  order  to  obtain  a  bilinear  relationship  between  any  two  solu- 
tions of  the  system  we  consider  the  real  adjoint  problem.  In  wave  propaga- 
tion problems  the  real  adjoint,  as  opposed  to  Hermitian  conjugate,  appears 
appropriate,  since  the  complex  conjugate  of  dispersion  function,  see  Eqs. 
(28)-(2S),  is  not  an  analytic  function.  The  equation  which  is  the  real 
adjoint  of  (41)  is 

W'  =  -M^  •  W,  (46) 

where   the   superscript  T  denotes  the  transpose.   If  T  is  a  solution  of  the 
original  equation,  and  W  is  a  solution  of  the  adjoint  equation  (46),  then 

Y*W  =  constant.  (47) 

After  examining  a  few  special  cases  of  the  adjoint  system  (46),  we 
found  that  there  is  a  simple  relation  for  our  system  that  connects  the  ad- 
joint problem  to  the  original  problem.  Specifically,  we  found  that  there  is 
a  nonsingular  constant  matrix  P  such  that  if  7  is  a  solution  of  (41)  [or  W 
is  a  solution  of  (46)]  then  P~^'Y  satisfies  (46)  [or  ?•¥  satisfies  (41)]. 
We  may  show  directly  the  existence  of  such  a  matrix  P.  We  assume  that  W 
satisfies  (46),  and  set 


-22- 


I  =  P  •  W.  (48) 

A  necessary  and  safficient  condition  for  P  to  satisfy  the  assumed  properties 
is  that 

M  •  P  +  P  •  M^  =  0.  (49) 

We  can  easily  construct  P  by  direct  calculation  by  making  use  of  the  proper- 
ty that  it  must  be  constant.  We  find  that  P  is  determined  uniquely  up  to  a 
multiplicative  constant,  and  that 


^N-k-l.k  "  P2N-k.k-l  =0-  f-  2  i  k  <  2N-2. 
^2N-1.2N  "  ~  ^2N.2N-1  "  ^   ^2N-2,1* 


To  make  the  solution  unique,  we  choose  P   ,«_,  =  1,  so  that 

^N-l-k,k  =  <-^^^'  ^  ik  <  2N-2. 
^2N-1,2N  "  ~  ^2N,2N-l  "  ~  ^' 


(50) 


(51) 


For   any   two  solutions  T  and  Z  of  (41),   the  bilinear  relation  can  then  be 
written  as 

N-1. 


Y  .  p-1  .  Z^:^^-l)^  <^j'^2N-l-j>  '^~' 


^^2N-1'^2N^  =  constant.     (52) 


where 


(Y  ,Z  )  =  Y  Z   -  Y  Z  . 
p   q      p  q     q  p 

We  now  obtain  the  desired  reciprocity  relations  by  selecting  the  solu- 
tions of  the  simple  scattering  problems  for  Y  and  Z  and  evaluating  them  by 
use  of  the  asymptotic  expansion  of  Y  and  Z  for  |s|  large  at  the  high-field 
and  low-field  sides.  Since  the  asymptotic  expansion  of  the  solutions  can  be 
decomposed  into  incoming  and  outgoing  fast  and  slow,  or  X  and  Q,  components. 
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it  is  sufficient  to  calcalate  the  bilinear  form  only  for  these  components. 
We  find  from  the  expansions  for  |s|  large  that  the  constant  value  of  the 
bilinear  relation  (52)  corresponding  to  the  pair  X  and  Q  is  zero,  namely 

X^  •  P  •  Q^  =  X^  •  P  •  Q^  =  0.  (53) 

The  superscripts  >  and  <  denote  incoming  and  outgoing  components,  respec- 
tively. For  the  fast  mode  the  contribution  of  the  sum  in  (52)  is  zero, 
giving 

X^  •  P  •  X^  =  -2in  lE^I^,  (54) 

whereas  for  the  slow  mode  only  the  sum  gives  nonzero  contribution  with  all 
summands  equal,  so  that 

Q>  .  P  .  Q<  =  -2i  f^i  IeJI^  (55) 

In  deriving  (55)  we  use  the  branch  of  (31)  corresponding  to  the  propagating 
Bernstein  mode  by  choosing  for  ts|  large 


.2N-2    -k   ,  .,k|   ,-k  .,,, 

A     =  P   =  (^i)  Ip  I   ,  (56) 

s  s 

and  we  note  that  the  first  component  of  Q  ,   namely  Q  ,   is  proportional   to 

Ip  1^^^,  cf.  (34). 

s 

We  now  apply  the  bilinear  form  to  the  solution  of  the  problems  of  X 
mode  incidence  from  high  and  low  field  sides.  These  solutions  for  |s|  large 
can  be  written  for  the  low  field  side  incidence  as: 

U^  =  af  X"^  +  b?  X^  LFS, 

I   >    I   >  <"> 

=  aj  X  +  c  Q^  HFS, 

and  for  the  high  field  side  incidence 
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U^^  =  bj^  X^  LFS, 

(  ^H'S 

II  >    II  <    II  > 
=  ag  X^  +  bg  X"-  +  c^^  q\  HFS. 

Here  we  use  the  notation  that  the  roman  numeral  superscripts  to  denote  the 
problems  I,  II,  III.  which  correspond  the  fast  mode  incident  from  HFS,  the 
slow  mode  incidence,  and  the  fast  mode  incident  from  the  LFS,  and  the  sab- 
scripts  L  and  H  to  indicate  the  low  and  high  field  side  components,  respec- 
tively.   By  evaluating  the  bilinear  form  in  both  sides  we  find  readily  that 

a.b,   =  a„b„  ,   which  states  that  the  transmission  coefficients  are  the  same 
L  L      an 

for  the  X  mode  incidences  from  the  either  side,  since 

In  order  to  obtain  the  reciprocity  relations  for  the  mode  conversion, 
we  consider  the  third  basic  scattering  problem  in  which  the  Bernstein  mode 
is  incident.   We  can  write  this  solution  for  |s|  large  as 

U^"  =  bf"  X<,  LFS. 

(60) 
III  _<  ^  III  ->  ^  .III  -<   ^„ 
=  a    X+c    Q+d    Q,  HFS. 

a 

When  we  apply  the  bilinear  form  to  problems  I  and  III,  and  evaluate  the 
constant  for  both  far  HFS  and  far  LFS,  we  obtain 

-.   I  .III    |_X,2     ,.   I  .III  Nzl  lcQ|2                      /ri\ 

-2i  a.  b,    ji  IE  I   =  -2i  c   d    t — i  |E^|  .                    (61) 

L   L        +                    I P  I  + 

s 

Multiplying  both  sides  by  their  complex  conjugates  and  rearranging  the  fac- 
tors we  easily  obtain  the  equality  of  the  mode  coupled  transmission  coeffi- 
cients for  the  fast  and  slow  mode  incidences.  This  establishes  the  second 
reciprocity  relation. 

Finally,  we  apply  the  bilinear  form  to  the  solutions  of  the  problems  II 
and   III.    In  this  case,   the  form  vanishes  at  the  far  low-field  side,   and 
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evalaation  at  the  far  high-field  side  leads  to 

,.   Ill  ^11    ,^i2   ,.   II  .III  N-1   ,_,Q,2  ,,,, 

0  =  2i  a^   ^e  ^  '  +'   ~   '^        Tp~I    +   ' 

s 

These  reciprocity  theorems  are  generalizations  of  the  results  of  Chin 
and  Man*  for  the  special  case  when  N=2 .  We  remind  the  reader  that  there  are 
two  more  exact  resnlts  previously  mentioned.  From  the  theorem  in  the  Appen- 
dix we  have  shown  for  N=2,  and  we  hypothesize  for  N>2,  that  (i)  the  trans- 
mission coefficient  for  the  X  mode  incidence  is  given  by  (40),  and  is  the 
same  as  that  given  by  geometrical  optics,  and  (ii)  there  is  no  reflected  X 
mode  for  a  high  field  side  incident  X  mode. 

IV.  FUNDAMENTAL  ION  RESONANCE 

We  obtain  field  equations  valid  for  this  case  by  substituting  (22)  into 
Maxwell's  equation,  to  get 

2 

-.2/6^  ,         1   2,„   ^  ,,  ^  1  _    3e'*  ^  .2     .2^1   2.2,_ 
^   ^^  ^i'l,0  -  2  %^%  ^  <1  *  2  ^i  -  6^  ^i"^  "2.0*^  ^  I  V  ^^- 

-irvS.i^^ii  =0- 

*^^-^-^i^t^^i'^^.2^^2'^^'>^l|  =0- 

where  we  have  again  used  the  charge  neutrality  n  =Z.n.  to  evaluate  the  elec- 

e   1  1 

tron  contribution  to  the  dielectric  tensor.   Equations  (63)  indicate  that  an 
appropriate   scaling  for  our  system  is 

5  =  e  , 

E^  =  0(e^E_).  E||  =  0(m^/mjE_). 
By  dropping  the  terms  of  order  m  /m, ,  and  eliminating  E  ,  we  obtain 
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[1  +  d^  +  d^A(s)d^]E  =  0.  (64) 


where  d  and  s  are  as  given  in  (24)  and  (25)  with  N=l,  and 

A(s)  =  T=-  '—r- 

2     ^     2   '-^ 
n  =2+X.  -nil. 
c         1     II 

We   identify  n    as  the  cold  plasma  refractive  index  for   the   magnetosonic 

(fast)   mode  evaluated  at  the  resonance.    We  may  generally  ignore  the   last 

term  in  the  definition  of  n   in  (64),   since,  as  we  show  by  numerical  compn- 

c 

tation   in  Sec. VI,   geometrical  optics  becomes  valid  whenever  nii  is   large 

enough  to  make  significant  to  make  a  significant  modification  to  n  . 

c 

We  may  extract  some  conclusions  from  a  simple  examination  of  (64). 
First,  although  we  have  carried  relativistic  modifications  of  the  Vlasov 
equation,  it  should  be  no  surprise  that  relativistic  effects  in  the  ions  are 
negligible.  Specifically,  the  relativistic  correction  included  in  the  fac- 
tor A.  defined  in  (4)  is  small  when 

n   ~  X.  >>  Ti 

c       1 

We  restrict  our  study  to  this  range,  so  that  a  functions  may  be  replaced 
by  the  standard  plasma  dispersion  functions.  Second,  the  Doppler  broadening 
effects  generated  by  the  term  containing  the  parallel  refractive  index  domi- 
nate spatial  gradients  in  the  boundary  layer  when 

s\   3  „2 

When  the  latter  condition  is  satisfied  we  expect  the  geometrical  optics  ap- 
proximation to  hold  within  the  resonance  layer.  We  verify  this  result  in 
our  numerical  investigation  which  we  present  in  Sec. VI.  Third,  since  the 
order   of  Eq.(64)  is  four  we  expect  mode  coupling  to  some  other  mode  to  oc- 
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cnr.  We  identify  this  mode  in  subsection  B  by  examining  the  uniform  medium 
dispersion  relation  near  the  resonance  layer.  Finally,  we  note  that  (64)  can 
be  solved  by  quadratures  for  the  case  of  purely  perpendicular  propagation  in 
the  nonrelativist ic  limit.  This  special  case  is  of  particular  importance 
since  it  also  corresponds  to  the  large  |s|  limit  of  Eq.(64).  We  present 
this  limit  in  the  following  subsection. 

A.  Asymptotic  Solotions 

For  |s|  large,  we  replace  A(s)  in  (64)  by 

A(s)  =  a/s. 

where 

a  =  (7<w'|>  -  <w?>^)n^6(Lo)/c)/(8X.).  (65) 

4  i      i    c  1 

For  a  Maxwellian  equilibrium  distribution  function  the  first  factor  in  the 
right-hand  side  of  (65)  is  1/2. 

Four  linearly  independent  solutions  of 

[1  +  d^  +  d^  (a/s)  d^]u  =  0 

are 

2  3         _. 

u  =  a  /   dp  -f-   eip[(s-a)p  +  ^  +  fr  log^lfl.  i^^) 

C.    p2^j  3    21    p+i 

where  the  contours  C,   i=l,2,3,4,   are  sketched  in  Fig. 1.   Contributions  to 

these   solutions  for  |s|  large  come  (i)  from  the  branch  points  p,  ^i,   and 

b 

2 
(ii)  from  the  stationary  points  p   =  (a-s)/a.   The  first  contribution  gener- 
ates  the  fast  branch  of  the  extraordinary  mode,   whereas  the  latter   corre- 
sponds  to   the  slow  branch  which  we  study  in  the  geometrical   optics   limit 

the  next  subsection.   We  omit  the  solution  corresponding  to  C  ,  which  grows 

4 

exponentially  in  the  low  field  side  and  thus  is  of  no  interest  in  the  solu- 
tions of  simple  scattering  problems.  For  the  remaining  three  solutions  we 
obtain  for  |s|  large: 
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Uq  =  e^*"''^  (2sinh^)(p^  +  p_)  +  (q^  +  q_) ,  HFS 


=  0,   LFS,  (67) 


5an/4 


a_  =  p,  0      ■*■  S'  ^^ 


=  p  e     ,  LFS, 


where 


P^  =  -  Iy  n^f7)exp[Ti(s-^a+|log2|s-2a|)], 

,  ,-   2  an 
I     |l/2  p  e  - 

an      _s .  ilL, j2i   i,    ., 

s      p^+  1 

The  basic  scattering  problems  in  this  limiting  case  can  be  constructed  as 
appropriate  linear  combinations  of  these  solutions.  First,  we  mast  investi- 
gate the  geometrical  optics  region  near  the  resonance  layer,  and  match  oar 
fast  and  slow  branches  represented  by  p  and  q  fanctions  given  above  to  those 
predicted  by  the  aniform  field  dispersion  relation.  We  carry  oat  this  anal- 
ysis in  the  next  section  in  which  we  also  determine  the  energy  flax  asso- 
ciated the  various  waves. 

B.  Geometrical  optics  region 

Consistent  with  the  formulation  given  in  subsection  A,  we  assume  that 
the  geometrical  optics  assumptions  hold  for  the  system  at  least  outside  of 
the  resonance  layer.  We  start  with  the  linearized  Vlasov  equation  with 
plane  waves  propagating  perpendicular  to  the  equilibrium  magnetic  field,  and 
calculate  the  dielectric  tensor  in  the  usual  manner.  We  employ  the  usual 
(+,-)  representation,  and  evaluate  the  tensor  elements  near  the  fundamental 
resonance  to  obtain 

%.  =  ^  -  ^i  ^  2^it°l,l  'h.2-%,0^^'-^^^'^- 
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where 

G.  ,  =  /d^w  W|  f^  J.  J^, 
i.j         i  3w|   i   j' 

J   are  the  Bessel  fanctions  with  argoments  miwiu/fl,   and  n  is  the  local  per- 
pendicular  index   of  refraction.    We  obtain  the  dispersion  relation   valid 

near  the  resonance  as 

2 

„2  ^  ,  ,  j^  ^.^  J  ^0,2  -  ^0,0^2.2  . 
1    l-a./Q   1     Gjj^Q 

If   we  keep  only  the  lowest  significant  order  terms  in  r\,      we  find  a   biqua- 
dratic dispersion  relation: 

4 

'^'^^  "^i  *  r:wi  ^i  i  "^'^H' '  ''^i'^-  ^''^ 

2  2 

We   readily   observe   that  if  we  replace  d   in  (66)  by  -(n/n  )    we   recover 

4 
(68).   For  j\   /(l-(o/Q)  small,   one  root  of  (68)  is  approximately  n  ,   corre- 
sponding to  a  fast  wave,  where 

4  =  %  "  ^  "^  ^i  ~  ^i-  <^9) 

4 

The  other  root  n  ,  again  for  n  /(l-u/fl)  small,  corresponds  to  the  slow  wave, 

2   2 
n  >>n.,  where 
s    f 

^2  ^  ^^j2  l-co/Q 

s      1    4 
e 

When  we  rewrite  (68)  in  the  more  convenient  form 

2    2  ^  a  n^ 

f   s   2 

•"f 


we  easily  see  that  the  slow  mode  propagates  only  on  the  high  field  side,  and 

2 


2   2 
coalesces  with  the  fast  mode  at  s=4a,  where  n  =2n,.   No  modes  propagate  from 


-30- 

this  valae  of  s  to  resonance.*  We  can  relate  the  refractive  indices  of  the 
slow  and  fast  modes  by 

2   s  2  ^   2  2  ,,^, 

%  =  r  °f  ~  -Ps  ""{•  (^o> 

Figure  2  illustrates  these  modes  in  the  (n.u/Q)  plane  for  a  fixed  value  of 
X.  and  T..  We  may  expect  that  some  portion  of  the  fast  wave  energy  incident 
from  the  LPS  should  tunnel  through  the  resonance  layer  to  the  HFS  and  par- 
tially mode  coupled  to  the  slow  mode.  Similarly,  HFS  incident  fast  or  slow 
modes  may  suffer  reflection  and  mode  conversion  in  addition  to  tunnelling  to 
the  LFS  as  a  fast  mode.  We  next  establish  the  relation  between  the  energy 
flux  associated  with  these  modes,  which  is  needed  for  the  scattering  prob- 
lems. 

The  component  of  the  energy  flux  along  the  direction  of  propagation, 
F  ,  for  our  problem  is 

F   =  !-  n   Im{E*dE  }. 

x   4jt  c    y  y 

2  i 

Since  E  =  0(e  E_) ,  we  may  replace  E  by  ~  E_,  and 

f'    n   Ie'I^        IeN^ 

F^    n,  lE^I^     ^  \eU^ 
X     f    -  - 

We  can  now  calculate  the  solutions  of  the  basic  wave  propagation  problems 
when  nil  =0. 

By  making  use  of  asymptotic  expansions  (67)  of  the  solutions  to  the 
wave  equation  (66)  we  can  readily  calculate  the  solutions  of  the  basic  scat- 
tering problems.  When  the  fast  mode  is  incident  from  the  HFS  the  solution  of 
the  problem  may  be  represented  by  u_.  We  find  that  the  transmission  coeffi- 
cient (or  the  fraction  of  the  transmitted  energy  flux)  for  the  fast  mode 
incidence  is 
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Tj  =  T  =  e    . 


(72) 


where  a  is  given  in  (65).  In  this  case,  there  is  no  reflected  fast  mode, 
and  the  energy  not  transmitted  is  fully  mode  converted  to  the  slow  mode. 
i.e.  R^  =  1-T.  To  obtain  these  results,  we  use  the  energy  flux  ratios  eval- 
uated in  (71). 

When  the  fast  mode  is  incident  from  the  LFS  the  solution  is  that  linear 
combination   of  u   and  u_  which  eliminates  p   on  the  HFS,   and  we  find   that 

the  transmission  coefficient  is  again  T  =  T,   the  reflection  coefficient  IL 

2 
equals  (1-T)  ,   and  the  mode  conversion  coefficient  is  T  =  T(l-T).   For  the 

M 

slow  mode  incidence  the  reflection  coefficient  is  R_  =  "i  ,  and  the  transmit- 
ted and  reflected  mode  conversion  coefficients  are  the  same  as  the  case  of 
fast  mode  incidence.  For  any  incident  wave  there  is  no  net  energy  absorp- 
tion for  nn  =  0.  For  non-zero  energy  absorption  we  must  assume  n n  ^  0,  and 
we  rely  on  numerical  computations  of  Sec.VT. 

C.  Reciprocity  relations 

We  write  (64)  in  the  form 


B'  =  M  .  u. 


(73) 


where  the  prime  denotes  derivative  with  respect  to  s,   the  components  of  the 
vector  n  are  u   =  E_,  u   =  u',  u   =  A(s)u'  +  u  ,  u   =  u'   and 


H  = 


0  1 
-1/A  0 
0  0 
-1    0 


0  0 

1/A  0 

0  1 

0  0 


One  can  easily  show  that  if  o  satisfies  (73),  then  v  =  P  •  n  satisfies  the 
real  adjoint  equation,  namely 


V'  =  -M 


▼, 


(74) 
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where 


[■  0-101 

P  =              I  10-10 

I  0      10     0 

L  -1    0    0    0 


We   hence   conclade  that  If  a  and  a  are  any  two  solatlons  of  (73)   then   the 
bilinear  form 

a  •  P  •  u  =  constant. 
The  reciprocity  relations  are  then  readily  obtained  by  direct  evaluation  of 
the  constant  at  the  far  LFS  and  far  HFS  for  any  pair  of  solutions  of  two 
distinct  basic  scattering  problems.  We  suppress  the  details  of  these  calcu- 
lation here.  We  find  that  the  transmission  coefficient  for  the  fast  mode 
incidence  is  independent  of  the  direction  of  incidence,  and  the  reflected 
(or  transmitted)  mode  coupling  coefficients  are  the  same  in  the  case  of  fast 
mode  converted  into  slow  mode  or  conversely.  In  proving  these  relations  we 
break  the  solutions  under  consideration  into  incoming  and  outgoing  fast  and 
slow  components,  and  utilize  their  asymptotic  forms  for  |s|  large,  which  we 
derive  in  the  next  section. 

The  theorem  of  the  Appendix  applied  in  this  case  shows  that  the  trans- 
mission coefficient  for  incident  fast  mode  is  given  exactly  by  (65),  which 
again  matches  the  geometrical  optics  result.  Also,  when  the  fast  mode  is 
incident  from  the  high  field  side  there  is  no  reflected  fast  mode. 

D.  Reformulation  of  the  general  case 

The  existence  of  solutions  growing  exponentially  for  |s|  large  in  the 
LFS  makes  it  necessary  to  use  some  special  numerical  integration  techniques. 
We  follow  closely  the  methods  described  in  detail  for  the  problem  of  second 
harmonic  electron  cyclotron  resonance.'-  We  first  transform  (73)  using  a 
Riccati  reduction  into  a  nonlinear  third  order  system  by  setting 


'i  =  ^+l/^•  i=i'2.3. 


(75) 
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This  transformation  essentially  reduces  the  exponential  growth  of  the  solu- 
tion into  algebraic  growth,  and  it  avoids  the  occurrence  of  very  large  num- 
bers.  We  find 


(76) 


'2  ^  '3  ~  '1*2' 

W^      =     -1     -     WgW^. 

We  integrate  (76)  from  far  LFS  (s  — >  -")  into  the  resonance  region  (s=0)  to 
obtain  one  solution  which  vanishes  asjnnptotically  in  far  LFS.  In  order  to 
obtain  two  other  solutions  which  behave  like  incoming  and  outgoing  fast  mode 
in  the  LFS  we  use  the  method  of  redaction  of  order  by  which  the  growing 
solutions  of  (74)  are  suppressed.  If  we  denote  by  y  the  solution  which  is 
exponentially  small  for  s  — >  -°>,  the  other  two  solutions  n  are  given  by  so- 
lutions z  of  the  reduced  system  by  the  relations 

u  =  z  +  Xy.  (77) 

z'  =  M  .  z  -  (z2/yj)y.  (78) 

where  we  have  set  z   =  0,  so  that  X  =  u./y^,  and  X'  =  z-/y-. 

For  s>0  all  solutions  of  (74)  represent  propagating  waves,  so  that  we 
may  return  to  (74).  We  extend  each  of  the  three  solutions  for  s<0  into  the 
domain  s>0,  and  for  convenience  we  choose  X  =  0  and  y  =  1  at  s=0.  We  inte- 
grate these  solutions  for  s  large  and  positive  (far  HFS) ,  and  we  decompose 
numerically  the  solutions  so  obtained  into  incoming  and  outgoing,  fast  and 
slow  waves.  The  solutions  of  basic  scattering  problems  can  then  be  obtained 
easily  by  taking  appropriate  linear  combinations  of  these  solutions.  Final- 
ly, by  making  use  of  the  energy  flux  theorems,  we  calculate  the  propagation 
characteristics   for   each  problem.    Numerical  results   are   presented   in 
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Sec.VI. 

V.  ION-ION  HYBRID  RESONANCE 

We  consider  a  plasma  containing  two  species  of  Ions  for  which  one  spe- 
cies is  at  fundamental  resonance  while  the  other  is  at  second  harmonic  reso- 
nance. To  be  more  explicit  we  refer  to  these  species  as  F  and  S  which  may 
represent  a  H-D  or  a  He-T  system.  Although  the  equations  for  the  species  F 
are  similar  to  those  given  in  the  previoas  section,  there  are  several  impor- 
tant differences  which  introduce  confusion.  We  define  the  boundary  layer 
variable  as 

5  =  M^  -  -^ ].  (79) 

6   2       0) 

or  equivalently 

Qt:(x/L) 
=  —  [1  -  — ] 

Note  that  (79)  is  to  be  contrasted  with  (2).  In  the  current  formulas  (18) 
and  (22).  (2)  is  used  to  define  5,  while  in  this  section  (79)  defines  5 
differently  for  each  species  at  resonance.  Thus,  we  must  be  careful  in  the 
evaluation  of  the  ion  currents.  We  introduce  the  parameters  r\„  and  r\„, 
where 


which  we  assume  to  be  on  the  same  order  of  magnitude,  and  we  assume 


"f.s  =  \.S  \   ^^   ^' 

where  n   is  the  appropriate  cold  plasma  refractive   index,   which,   together 
c 

with  the  stretching  factor  8,   we  specify  later.   We  find  the  ion  current  to 
be   the   sum  of  (18)  and  (22).    We  readily  observe  that  our  boundary   layer 


-35- 


expansion  is  consistent  provided  that 

8  =  e„, 

,  (81) 

Xp/Xg  =  (Zp/Zg)''(nj,/ng)(mp/mg)  =  2(Zg/Zp)  (n^/ng)  =  0(8). 

We  add  the  electron  current  obtained  from  the  charge  neutrality,   drop  terms 

of  order  m  /m.,  and  eliminate  E  ,  to  obtain 
e   1  - 

{1    +   d^    +    (1    +   H   +   d^)[dA(s)d   -   B(s)]}E^   =   0,  (82) 


where 


2  2 

\i  =    (e   _  -   nji  )/  (e^,^  ~  °||  ^  ' 

\  =  (e__  -  4 )(%+  -  4  )/(e„  -  4 ). 

d  =  d/ds. 


s   =  -2£8n   Ixo/c. 
^      c 

^^^^    =^s"l,o/^^x-4>' 


B(s)    =^^. 


e      -n 

XX 


II 


According  to  (81)  the  density  ratio  n^/n„=0(8)  and  this  might  seem  to 
be  an  uninteresting  parameter  range.  However,  the  numerical  calculations 
given  in  Sec. VI  show  that  a  maximum  in  the  wave  energy  absorption  occurs 
with  a  rather  small  addition  of  minority  species  on  the  order  of  1%  of  the 
electron  number  density.  Such  a  small  minority  fraction  is  consistent  with 
our  ordering.  We  note  that  if  the  ordering  (81)  is  not  satisfied,  within 
the  layer  E  must  get  much  smaller  in  magnitude  in  order  to  keep  the  first 
term  in  the  right-hand  side  of  the  first  equation  in  (22)  balanced.   We  find 
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that   in   this  case  a  second  boundary  layer  exits  with  5=8    within  which 

2  2 

(64)  is  satisfied,  with  n  =  2  -  L  +  ^^S^^  ~  ^°II  *  ^°^  *^^  practical  pur- 
poses, however,  the  inner  boundary  layer  is  almost  transparent,  since  the 
transmission  coefficient  determined  by  (72)  is  very  nearly  one. 

Ve  have  not  found  explicit  integral  representations  of  solutions  of  the 
asymptotic  form  of  (82)  for  |s|  large.  We  may  however  extract  necessary 
informations  about  the  asymptotic  behavior  of  solutions  of  (82)  by  making 
use  of  standard  techniques.   We  write  (82)  in  the  form 


n'  =  M  •  u. 


(83) 


where  the  prime  denotes  derivative  with  respect  to  s,   the  components  of  the 

vector   n   are   u,    =   E    ,    u.    =   A(s)u',    u,    =   u,    +   u'    -   B(s)u, ,    u.    =   u ' ,    and 
1+2  1312  143 


t-1 

0 

1 

0 

0 

0 

0 

1 

n 

0 

-l-^ 

0 

(84) 


If  we  set  w,  =  u.  ,/a, .   j=l,2,3,  we  find  the  Riccati  reduced  system  for  our 


J 


problem  may  be  written  as 


(85) 


w^  =  w^  -  1  +  B(s)  -  w^/A(s), 

w'  =  w   -  w  w  /A(s) , 

w'  =  (1  -  (l+(i)w  -  w  w  /A(s). 
We   construct  asymptotic  expansions  of  solutions  of  (85)  in  terms  of   asymp- 
totic series  in  reciprocal  powers  of  s.   For  the  X  mode,  we  find 


w^  =  TiA  [1  -  ^(A+B)], 
w^  =  1  -  A  -  B. 
w^  =  Ti  [1  -  (1  -  ^)(A+B)], 
where  A  and  B  denote  a/s  and  b/s,  respectively,  and 


(86) 
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X„   6n   L<d/c 
b        c 

*   =  -  ~2 

"xx-^ll 


(87) 
X^,  n  L(i)/c 
b  =  -  -^^ 


2    , 


For    the   Bernstein  mode 


w,    =   .iA^/^d    *  JA  -  iB)    .  ^. 

3/2  ^«»> 

w     =  -^lMl   +  A  +  B)    ±   ijiA'^"'/s, 

w^    =  ±inA^^^[l    +    (1    +  ^)A   +  i-B]    -    n^    . 

Since   u'/u  e  w  /A  is  proportional  to  the  derivative  of  the  eikonal  function 

for  |s|  large,  we  readily  find  for  the  X  mode 

u  =  e'''^^'  "  (n/2)(a+b)log  s]^  ^^^^ 

For  the  Bernstein  mode,  we  obtain 

u^  =  s^/*  e^i^"'''(j^^^»-^\  (90) 

The  theorem  of  the  Appendix  applied  to  this  case  shows  that  the  trans- 
mission coefficient  for  the  X  mode  is 

T^  =  e-'^^^'^*^^'".  (91) 

Furthermore  the  reflection  coefficient  is  zero  for  the  X  mode  in  high  field 
side  incidence.  We  next  match  these  resolts  to  the  geometrical  optics  re- 
sults evaluated  near  the  resonance  layer,  and  establish  the  energy  flux  re- 
lations. 

A.  Geometrical  optics  results 

For  the  present  problem,  we  can  approximate  the  dispersion  relation  for 
the  Bernstein  mode  as  e  =0,  which,  when  evaluated  near  the  resonance, 
yields  the  relation  between  the  refractive  index  of  the  quasi-electrostatic 
(Bernstein)  mode  and  n  : 
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2 

!2 


A  "^  5X 


h     . 


(92) 


Eqaation  (92)  also  is  a  direct  consequence  of  (90). 

When  we  evaluate  the  component  of  the  energy  flax  along  the  direction 
of  propagation  for  the  Bernstein  mode  we  find  that  the  contribntion  from  the 
species  F  is  negligible  near  the  resonance  layer,  so  that 


,Q 


.Q,2 


F^  lE^I^ 


(93) 


We   use    (93)    to   calculate   the  mode   coupling   coefficients   associated  with   the 
scattering   problems. 

Before  we  present  numerical  results,  we  give  some  analytical  results 
directly  obtained  from  (82).  First  we  prove  that  the  reciprocity  relations 
are    still    true    for    the    ion-ion   hybrid   case   discussed   here. 


B.    Reciprocity    relations    for    the     ion-ion    hybrid    resonance 

The      proof   of    the    reciprocity   relations   can  be    established    in   a     manner 
analogous    to   the   previous   case.      We   find   that    if 


P  = 


0  n     0     0 

-H  0      0      0 

0  0     0      1 

0  0-10 


then  if  a  satisfies  (83)  then  v=P-o  satisfies  the  real  adjoint  equation 

T 
V'  =  -M  -v. 

Hence,  if  u  and  a  any  two  solutions  of  (83),  then  the  bilinear  form 

U'P-u  =  constant. 

By   evaluating  the  constant  at  far  LPS  and  far  HFS  for  the  solutions  of   the 

basic   scattering  problems,   we  again  derive  the  same  reciprocity   relations 
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as  for  the  other  resonances. 

VI.  NUMERICAL  RESULTS 

A.  Absorption  at  the  third  cyclotron  harmonic  resonance 

The  results  obtained  in  Sec.III  indicate  that,  because  b  in  the  expres- 
sion (40)  is  typically  large,  when  the  propagation  is  purely  perpendicular 
to  the  equilibrium  magnetic  fieldfast  waves  incident  the  Nth  ion  cyclotron 
harmonic  resonance  layer  are  mainly  reflected  for  the  low  field  side  inci- 
dence, and  converted  to  reflected  ion  Bernstein  waves  for  the  high  field 
side  incidence.  When  the  propagation  is  not  purely  perpendicular  exact  so- 
lutions are  not  available,  and  we  turn  to  numerical  computation.  In  this 
section  we  study  wave  propagation  near  third  harmonic  resonance.  We  find 
that  for  the  low  field  incident  fast  mode,  when  nii  is  approximately  5  or 
larger,  reflection  and  mode  conversion  of  this  mode  become  small,  so  that 
all  incoming  wave  energy  which  is  not  transmitted  to  the  high  field  side  is 
absorbed  by  the  plasma.  The  rate  of  decrease  in  the  reflection  and  mode 
conversion  coefficients  as  a  function  of  n •■  increases,  and  the  correspond- 
ing increase  in  the  absorption  coefficient,  is  largely  independent  of  the 
plasma  temperature  or  density.  However,  this  rate  increases  with  the  charac- 
teristic length.  We  thus  conclude  that  these  waves  may  provide  an  efficient 
way  of  heating  a  plasma,  particularly  in  a  large,  moderately  hot  tokamak. 

For  waves  incident  from  the  high  field  side,  as  a  consequence  of  the 
much  slower  rates  of  decrease  of  the  reflected  mode  conversion  coefficient 
with  increasing  parallel  refractive  index,  the  absorption  rates  are  much 
smaller.  When  Um  becomes  sufficiently  large,  the  absorption  is  further 
reduced,  since  the  transmission  coefficient  becomes  larger  due  to  geometri- 
cal effects.  The  absorption  coefficients  for  the  high  and  low  field  inci- 
dent fast  mode  eventually  equal  one  minus  transmission  coefficient,  and  geo- 
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metrical  optics  becomes  valid.  The  valae  of  Qii  in  this  region  is  quite  high 
so  large  that  accessibility  of  waves  may  be  poor.  For  the  case  of  Bernstein 
mode  incidence  we  find  that  the  absorption  is  poor  when  n,,  is  small,  which 
is  mainly  caused  by  the  large  reflected  mode  conversion.  and  the  absorption 
sharply  approaches  one  when  Um  is  sufficiently  large. 

In  our  numerical  study,  we  utilize  a  two-fold  reduction  of  order  tech- 
nique to  eliminate  the  two  linearly  independent  solutions  which  grow  on  the 
low  field  side.  By  integrating  towards  the  resonance  we  calculate  numerical- 
ly two  solutions  which  are  zero  asymptotically  .  In  addition  we  also  cal- 
culate two  (incoming  and  outgoing)  fast  mode  components.  On  the  high  field 
side  there  is  only  one  solution  which  grows,  and  we  eliminate  it  by  the  same 
method  given  in  our  treatment  of  second  harmonic  resonance. ^  By  inte- 
grating again  towards  the  resonance,  we  calculate  one  solution  which  is  as- 
ymptotically zero  and  four  solutions  which  are  propagating  fast  and  slow 
modes.  We  match  the  solutions  computed  on  the  low  field  side  to  the  solu- 
tions computed  on  the  high  field  side  at  the  resonance,  and  thus  obtain 
solutions  valid  within  the  entire  boundary  layer,  giving  the  solutions  of 
the  basic  scattering  problems.  We  also  use  the  reciprocity  relations  ob- 
tained in  Sec. Ill  to  improve  the  numerical  accuracy  of  the  results  which  we 
present  next. 

In  Fig. 3  we  present  the  wave  propagation  characteristics  in  a  hydrogen 


p 


plasma   for  various  values  of  X  =  ((u^/w)   as  a  function  of  nn  ,   for  T  =  10 


keV,  and  Ixo/c  =  1,  where  L  is  the  characteristic  length.  The  results  just 
stated  are  clearly  evident.  We  see  that  as  nii  increases,  the  transmitted 
mode  conversion  coefficient,  T„,  and  reflection  coefficient  for  the  low 
field  incident  fast  mode,  R^ ,  sharply  approach  zero.  A^ ,  the  absorption 
coefficient  for  the  low  field  incident  fast  mode,  then  equals  1  -  T  ,  as 
expected  by  geometrical  optics.    The  rate  of  decline  in  the  reflected   mode 
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coupling   coefficient,   R^,    as  a  function  of  nn  is  much  slower  than  T^  and 

IL  ,   and  it  depends  on  the  value  of  X.    As  a  result  of  this  the   absorption 

coefficient   for  the  high  field  incident  fast  mode  is  more  sharply  peaked  at 

a  much  lower  value  that  that  of  the  low  field  incidence.    The   absorption 

coefficient  A   for  the  Bernstein  mode  incidence  increases  monotonically  with 
o 

nn  .   We  find  similar  behavior  for  other  values  of  T. 

Figure  4  depicts  the  absorption  contours  for  the  low  and  high  field 
incident  fast  mode  in  the  (X.nii  )  plane.  We  also  present  similar  results  in 
the  (T»°||  )  plane  in  Fig, 5.  We  study  the  effect  of  the  machine  size  on 
absorption  in  Fig. 6.  We  see  that  as  L<i>/c  increases  the  mode  coupling  ap- 
proaches zero  more  rapidly  with  nn  ,  which  results  in  higher  absorption 
rates. 

B.  Wave  propagation  at  the  fundamental  ion  resonance 

In  Sec. IV  we  have  shown  that  the  boundary  layer  method  can  be  applied 
successfully  to  case  of  fundamental  resonance,  where  within  a  narrow  reso- 
nance region  the  fast  and  slow  modes  bifurcate  and  geometrical  optics  fails. 
We  obtained  an  analytical  expression  for  the  transmission  coefficient  for 
the  case  of  fast  mode  incidence.  Reflection  and  mode  conversion  coeffi- 
cients are  also  given  for  the  case  of  purely  perpendicular  propagation.  In 
the  present  subsection  we  numerically  study  the  case  when  nn  is  not  zero  for 
an  hydrogen  plasma.  We  show  that  reasonable  absorption  rates  of  the  inci- 
dent fast  mode  energy  is  possible  if  the  ion  temperature  is  greater  than  10 
Key,  and  density  is  sufficiently  high  (DIOOO),  when  n  •■  >  2. 

Since  the  field  equation  is  similar  to  that  studied  for  the  second 
electron  resonance  which  we  have  studied  in  Refs.1-2,  we  employ  the  same 
numerical  method  in  solving  (73).  We  find  that  for  T  <^  5  keV  the  incident 
fast  magnetosonic  mode  largely  tunnels  through  the  resonance  layer,   causing 
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little  absorption,  even  for  densities  so  high  that  X  =  10000.  For  T  =  15 
keV  respectable  absorption  rates  may  be  achieved  at  larger  values  of  1  and 
moderate  valaes  of  nii  .  We  present  these  results  in  Fig. 7  where  the  absorp- 
tion contours  in  the  (I,n||  )  are  plotted  for  T  =  10  kev.  Figure  8  depicts  a 
similar  plot  for  T  =  15  keV.  Geometrical  optics  is  valid  when  the  absorp- 
tion rates  for  the  HFS  and  LFS  incidences  approximately  equal  each  other, 
namely  A^  ~  A. . 

We  find  quite  respectable  absorption  rates  for  the  incident  slow  mode 
even  when  both  density  and  parallel  refractive  index  are  small.  We  thus 
conclude  that  the  latter  may  be  considered  as  a  potential  heating  scheme  for 
tokamak  plasmas  if  these  waves  are  accessible  at  the  resonance  layer. 

C.  Numerical  study  for  the  ion-ion  hybrid  resonance 

We  have  shown  in  Sec.V  that  when  two  ion  species  are  present  in  the 
plasma,  one  at  the  fundamental,  and  the  other  at  the  second  harmonic  reso- 
nance at  the  same  resonance  surface.  the  energy  transmission  coefficient 
given  by  (92)  decreases.  Our  formulation  is  also  valid  in  the  limit  when 
the  minority  ion  density  approaches  zero,  so  that  the  numerical  results  pre- 
sented in  this  section  contains  as  a  special  case  the  results  with  only  one 
ion  species  resonant  at  the  second  harmonic.  Although  we  have  analytic  re- 
sults for  the  latter  case  when  nii  is  zero,  we  still  have  to  carry  out  a 
numerical  study  for  the  two-ion  problem. 

When  the  minority  density  becomes  sufficiently  large,  the  resonant  term 
in  J  is  balanced  by  E  ,  whose  magnitude  becomes  small  as  waves  get  close  to 
the  resonance,  and  waves  transfer  into  the  magnetosonic  mode.  As  we  have 
indicated  in  Sec.V,  a  second  boundary  layer  is  formed  in  this  case  which  is 
much  closer  to  the  exact  resonance  surface.  The  result  of  the  previous  sub- 
section  shows   that  incident  fast  waves  very  nearly  all  tunnel  through   the 
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resonance,  unless  X^  is  extremely  large.  We  thus  may  ignore  the  effect  of 
the  inner  layer  when  tokamak-like  plasmas  are  being  studied.  For  every 
tokamak  configuration  that  we  have  studied  the  maximum  increase  in  the  ener- 
gy absorption  rate  caused  by  the  minority  species  occurs  when  the  minority 
ion  density  is  somewhat  less  than  one  percent  of  the  electron  density.  The 
Bernstein  mode  is  not  affected  by  the  inner  boundary  layer.  Hence,  the  ef- 
fect of  large  concentrations  of  fundamental  species  need  not  be  considered. 
We  present  our  results  in  Figs. 9-11,  where  the  absorption  contours  are 
plotted  in  the  plane  of  minority  density  and  the  parallel  wavenumber.  We 
have  selected  three  examples  corresponding  to  typical  PLT,  JET,  and  CIT  pa- 
rameter  ranges.   Figure  9  depicts  the  case  of  a  (H-D)  plasma  with  B   =  3  T, 

13    -3 
n  =  3  X  10    cm   .   T„  =  T^  =  2  keV  and  L  =  120  cm.    We  see  that  although 
e  n    U 

the   fast  mode  incidence  yields  rather  poor  absorption  rates,   the   incident 

Bernstein  mode  can  be  used  to  heat  the  plasma  quite  efficiently.   In  Fig. 10, 

13 
we   give  the  results  for  a  (H-D)  plasma  with  B   =  3.45  T,   n   =  3.5  X   10 

\)  c 

cm   ,  T  =  T„  =  5  keV  and  L  =  300  cm.   We  find  good  absorption  rates  for  the 
n     D 

low   field   side  incident  fast  mode  even  for  quite  small  values  of   parallel 

wavenumber  when  n   is  about  .Oln  .   The  absorption  of  the  incident  Bernstein 
n  e 

mode  is  appreciable  at  moderate  values  of  kii  .    Finally,   in  Fig. 11  we  pre- 

3 
sent  the  absorption  level  lines  for  a  (  He-T)  plasma  with  B  =  10  T,   n  =5 

u  e 

X  10^^  cm~^,   T„   =  T„  =  10  keV  and  L  =  122  cm.    We  find  that  if  the  waves 
ne     1 

are  incident  from  the  high  field  side  they  suffer  a  nearly  full  reflected 
mode  conversion,  since  R^  ~  1.  We  find  excellent  absorption  rates  for  fast 
waves  incident  from  the  low  field  side  when  the  minority  density  is  0.5%  of 
the  electron  density.  More  energy  is  deposited  to  the  majority  ions  at 
larger  values  of  parallel  wavenumber.  We  also  give  the  deposition  profiles 
for  the  ion  species  within  the  resonance  layer  in  Fig. 12  for  a  typical   case 
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of  kii  =   2m   and  n—     =   0.005  n  .    We  see  that  the  energy  deposition   is 
qaite  localized  for  this  case. 

The  higher  absorption  rates  that  can  be  obtained  for  larger  valnes  of 
the  parallel  wavennmber  are  due  to  Doppler  effects,  rather  than  resonance 
effects,  hence  this  domain  should  be  avoided  when  accessibility  is  a  prob- 
lem. On  the  other  hand,  if  kn  is  too  small,  we  find  that  most  of  the  energy 
is  deposited  to  the  fundamental  species. 

VII.  DISCUSSION  AND  CONCLUSIONS 

This  paper  contains  analytical  and  computational  results  dealing  with 
wave  propagation  across  the  Nth  ion  cyclotron  resonance  layer  where  the  geo- 
metrical optics  approximation  is  not  valid.  We  have  used  a  boundary  layer 
expansion,  and  when  it  is  valid  we  can  solve  the  linearized  Vlasov  equation 
order  by  order  within  the  resonance  layer.  We  find  that  the  Nth  Fourier 
component  (N>^2)  in  gyrophase  angle  of  the  perturbed  distribution  function 
becomes  large  in  magnitude  within  the  resonance  layer,  so  that  there  is  a 
zeroth  order  contribution  to  the  current  density  from  the  terms  involving 
the  spatial  derivatives.  When  N=l ,  the  current  contribution  is  known  to  be 
large.  An  explicit  evaluation  of  the  resonance  layer  current  gives  rise  to 
field  equations  valid  within  the  layer  for  the  cases  of  pure  Nth  harmonic 
resonance,  NH ,  and  for  the  case  of  ion-ion  hybrid  resonance.  In  most  cases 
with  nil  =  0  we  may  represent  the  solutions  by  contour  integrals,  in  all 
cases  we  may  determine  the  transmission  coefficient  for  the  incident  fast 
wave  exactly,  cf.  F,qs.(40),  (72),  and  (92). 

For  pure  Nth  harmonic  resonance,  N>^2,  and  for  a  given  value  of  the 
effective  expansion  parameter  8K(0),  the  optical  thickness  (or  the  transmis- 
sion coefficient  T)  takes  a  minimum  value  (correspondingly  a  maximum  for  T) 

when  X.  satisfies  the  cubic  equation 
1 
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nV  +  2N^(3-3N-N^-N^)X^  +  (8N-9)N^(N^-1)X  -  2(N-2)  =  0. 

For  N  =  2  this  equation  has  three  real  roots,   0,  0.18,  7.32,  and  for  larger 

2 
values  of  N  the  real  root  is  approximately  2(N  +  1  +  3/N  -  3/N  ).   From  this 

result  we  easily  infer  that  when  N  is  greater  than  2,  for  almost  every  case 
of  physical  interest  (with  the  possible  exception  of  some  astrophysical  ap- 
plications where  L  may  be  very  large),  T  is  nearly  zero.  For  example,  for  N 
=  3  we  find  T  ~  exp(-429 .4/ |K(0) I )  =  exp(-429.4e^k^L) .  where  k^  =  Nu  /c  is 
the  cold  plasma  wavenumber  for  the  X  mode.  Therefore,  the  fast  mode  inci- 
dent from  the  high  field  side  is  almost  fully  mode  converted,  whereas  there 
is  an  almost  full  reflection  of  this  mode  when  the  incidence  is  from  the  low 
field  side.  This  result  clearly  indicates  the  breakdown  of  the  geometrical 
optics  expansion  near  the  resonance  layer. 

In  Secs.III-V  and  the  Appendix  we  have  presented  analytical  proofs  of 
various  properties  of  the  scattering  problems.  We  have  shown  that  reciproc- 
ity relations  hold  within  the  limits  of  our  expansions.  These  relations 
state  that  the  transmission  coefficient  is  independent  of  the  direction  of 
the  incidence  even  when  the  relativistic  corrections  and  the  effects  of  the 
nonzero  nii  are  included.  Further,  the  mode  converted  transmission  coeffi- 
cients for  the  fast  and  slow  mode  incidences  are  equal.  Similarly,  the  mode 
converted  reflection  reflection  coefficients  for  the  fast  and  slow  mode  in- 
cidence are  equal.  These  properties  are  very  useful  in  simplifying  numerical 
matching,  and  they  can  also  be  used  to  check  the  accuracy  of  the  numerical 
results.  We  have  also  proven  a  more  difficult  proposition,  which  gives  ex- 
plicit values  for  the  transmission  coefficient  for  the  incident  fast  waves, 
and  shows  that  there  is  no  reflection  for  the  fast  mode  incidence  from  the 
high  field  side.  These  proofs  are  given  only  for  the  special  cases  of  fun- 
damental, second  harmonic,  and  ion-ion  hybrid  resonance.  For  higher  ion 
harmonics   the  proof  becomes  too  complicated  to  be  considered  here   although 
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we  believe  that  the  proposition  is  still  trae. 

Section  VT  contains  the  results  of  nnmerical  studies.  When  nn  is  suf- 
ficientlj  large  we  recover  geometrical  optics  since  the  coefficients  for 
wave  reflection  and  moda  conversion  to  Bernstein  mode  become  small.  In  this 
limit  the  energy  absorption  becomes  approximately  one  minus  transmission.  We 
first  considered  the  special  case  of  N=3 .  On  the  basis  of  this  study  we 
conclude  that  waves  propagating  with  frequencies  equal  to  the  third  ion  cy- 
clotron harmonic  provide  an  efficient  means  for  heating  the  plasma,  particu- 
larly so  for  the  case  of  the  low  field  incident  fast  mode  with  a  small  nii  . 
We  have  also  shown  that  even  smaller  values  of  nn  are  needed  for  larger 
tokamaks.  Hence  one  may  achieve  better  accessibility  of  incident  waves  in 
this  case.  The  absorption  of  the  HFS  incident  fast  mode  and/or  the  Bernstein 
mode  incidence  require  considerably  higher  values  of  Un  in  order  to  achieve 
reasonable  absorption,  thus  they  are  not  as  promising  as  the  LFS  incidence. 
We  also  find  that  the  absorption  is  more  efficient  for  larger  ion  tempera- 
tures. 

The  case  of  magnetosonic  waves  incident  on  the  fundamental  resonance 
layer  when  nn  is  not  zero  is  also  studied  numerically.  We  have  found  that 
reasonable  absorption  rates  of  the  incident  fast  mode  energy  is  possible  if 
the  ion  temperature  is  greater  than  10  Kev,  and  density  is  sufficiently  high 
(X>>1000),  when  nii  >  2.  We  find,  however,  that  in  the  case  of  the  slow  mode 
incidence  reasonable  absorption  rates  may  be  achieved  for  the  parameter 
range  comparable  to  tokamaks  which  exists  today. 

In  the  case  of  two  ion  hybrid  resonance,  we  find  that  very  promising 
improvement  in  absorption  is  obtained,  particularly  when  the  minority  densi- 
ty is  approximately  one  percent  of  the  electron  density.  At  this  value 
there   is   a  strong  peak  in  the  absorption  rate  for   LFS  incidence  when   the 
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parallel  refractive  index  is  safficiently  small.  We  have  also  calculated 
energy  deposition  profiles  for  each  species.  Vhen  the  plasma  parameters  are 
chosen  appropriately,  we  find  that  the  plasma  heating  by  this  method  is 
quite  favorable. 
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APPENDIX 

In  this  appendix  we  prove  a  property  of  asymptotic  expansions  of  solu- 
tions of  differential  equations  appropriate  to  the  cases  of  fundamental, 
hybrid,  and  second  harmonic  resonance.  From  this  result  it  is  possible  to 
obtain  an  explicit  formula  for  the  transmission  coefficient  for  an  incident 
fast  wave,  and  to  show  that  for  a  fast  wave  incident  from  the  high  field 
side  no  fast  wave  is  reflected.  An  extension  of  the  proof  could  no  doubt  be 
given  for  any  higher  harmonic  resonance,  but  we  do  not  explore  this  possi- 
bility. 

We  consider  a  system  of  four,  linear,  first  order  differential  equa- 
tions, which  we  write  in  the  form 

y'  =  A(s)  y,  (Al) 

where  A(s)  is  a  four  by  four  matrix  and  y(s)  is  a  four  vector.  We  assume 
that  A(s)  is  analytic  in  Im  s  <  0,  and  continuous  in  Im  s  <.  0.  All  diagonal 
elements  of  A(s)  vanish,  and  more  generally  Tr  A(s)  =  0.  Additionally,  A(s) 
possesses  an  asymptotic  expansion  for  |s|  large  valid  in  Im  s  <,  0 

OD 

A(s)  =  s  5  a  /s°,  (A2) 

n=0 
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so  that  for  any  integer  M  there  is  a  constant  8L.  such  that 


|A(s)  -  s  >  a  /s^^l  <  K„|s|  "  (A3) 


n=0  '^ 

All   the   matrices  a   are  real  and  the  diagonal  elements  of  all   a   vanish. 

n  n 

All   solntions  of  (Al)  are  analytic  in  Im  s  <  0  and  continnons  in  Im  s  <^  0. 

In  addition,  for  any  four  vector  solutions  of  (Al),  say  ,  z, ,  z. ,  z, ,  and  z. 

I    2    J        4 

the  matrix  M  =  (z,  z.  z,  z.)  has  a  constant  determinant. 
12   3   4 

For  equations  of  the  form  (Al),  where  asymptotic  expansions  of  A(s)  of 
the  form  (A2)  are  assumed,  it  is  possible  to  obtain  formal  asymptotic  expan- 
sions of  solntions  which  are  of  the  form  of  the  product  of  a  power  of  s,  the 
exponential  of  a  polynomial  in  a  fractional  power  of  s,  and  a  reciprocal 
power  series  in  a  fractional  power  of  s.  Such  asjrmptotic  expansions  are 
developed  in  Coddington  and  Levinson.»  It  is  also  shown  there  that  on  any 
ray,  arg(s)  fixed,  there  are  solutions  of  the  original  system  which  are  as- 
ymptotic to  these  series  which  have  been  formally  obtained.  The  asymptotic 
series  are  uniquely  determined  by  the  coefficient  matrices,  and  conversely, 
if  the  formal  series  are  given  to  some  order,  then  the  coefficient  matrices 
are  also  known  to  some  order. 

Appropriate  to  the  systems  of  differential  equations  of  interest  to  us, 
we  assume  that  there  are  four  formal  asymptotic  series  of  the  form 


y  =  exp{is  +  i(b/2)log  s}   ^  ^J^^' 


(A4) 


n=0 


y^  =  exp{-is  -  i(b/2)log  s}  ^  0*/s°.  (A5) 


n=0 


1/4     ,      3/2  ^  .   1/2,  f    , 
y.  =   s     exp(cs     +ds    ))y/s 


(n/2) 


(A6) 


1/4     ,    3/2    .   1/2,  T  .  ,  (n/2)  ..-v 

y.  =  s     exp(-cs     -  d  s    )  >  8  /s      ,  (A7; 

n=0 

where  b,   c,  d.  are  real,  and  b  and  c  are  positive.  In  addition  p  ,  y  ,  and 
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8-   are   nonzero.    For  any  formal  series  of  the  form  (A2)  or   (A4)-(A7)   we 

introduce   the  notation  of  superscript  P  to  denote  the  series   truncated   to 

_p 
terms  of  order  s    times  any  possible  exponential  factor,  thus 

A   =  s  >  a  /s  , 

Jo   ° 


y   =  exp{is  +  i(b/2)log  s}  ^   P^/*  » 


n=0 


P     1/4     ,   3/2  ^  J   1/2,  ^   , 
y   =   s     exp(cs     +  d  s    )  )  T  /s 

n=0 


(n/2) 


and   so  on.    On  any  ray  arg(s)  fixed  there  are  four  solutions  of  (Al)  which 

P    P    P    P 

are  asymptotic  to  y  ,   y. ,   y, ,  y..   For  the  exact  solutions  it  follows  that 

the  determinant  of  the  matrix  of  the  four  solutions  is  constant.    Thus,   we 
conclude  that  there  is  a  constant  K' 


K^  <  ldet(y^  y^  yj  yj) I  <  1/K^ 


(A8) 


uniformly  in  Im  s  <_  0,  |s|  >_  1.  Here,  as  elsewhere  in  this  proof,  we  must 
exclude  a  neighborhood  of  the  origin,  and  consider  domains  lying  in  |s|  >.  1 
only. 

We  prove  that  there  is  a  solution  yp(s)  of  (Al)  defined  in  the  half 
plane  Im  s  <^  0  such  that  for  -5n/6  <.  arg  s  <_  0  y^  has  the  same  asymptotic 
expansion  as  y  ,  or  for  any  P  there  is  a  constant  Kp 


lyp  -  y^l  <  Kp  |exp{-is  -  i(b/2)log  sll/lsl^""^ 


(A9) 


while   for  -n  ^  arg  s  <.  -Sn/6   y_  has  the  same  asymptotic  expansion  as  y_   + 
Xpy,  or  for  any  P  there  are  constants  Xp  and  Cp  such  that 


lyp  -  y2  -  Vs'  ^  S  l"Pt-i«  -  i(b/2)log  s)|/|s|^'^V 


(AlO) 
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It  then  follows  that  j„  represents  a  fast  wave  incident  from  the  high  field 
side  (arg  s  =  0) ,  that  no  fast  wave  reflected,  and  that  the  energy  flax 
transmission  coefficient  is  exp(-bn)  .  The  mode  converted  amplitude  \  can- 
not be  determined  explicitly. 

The  proof  proceeds  in  fonr  steps.    First,  we  construct  an  equation  ap- 

P    P    P    P 

proximately  equal  to  (Al)  with  solutions  y, ,   y. ,   y. ,   y.  valid  in  the  full 

12     3     4 

domain   Im   s   <.  0.    We  solve  (Al)  by  perturbation   about   the   approximate 

equation   and  we  thus  formulate  an  integral  equation  whose   solutions   solve 

(Al)   and  which  are  also  close  to  our  approximate   equation.    Secondly,   we 

consider  the  integral  equation  in  the  sector  -Sn/6  <.  arg  s  ^  0  and  show  that 

it   has   solution,   called   y„,   which  has  as  its  assmptotic   expansion   y  . 

Third,   we  examine  the  sector  -n  <.  arg  s  <^  -2n/3  -  5,   0  <_  S  <.  n/6,   and   we 

show   that   there   are   four  solutions  of  the  integral   equation   which   are 

asymptotic  to  y  ,   y  ,   y  ,  and  y . ,  respectively.   The  fourth,  and  the  last, 
12     3        4 

element  of  the  proof  is  trivial  given  the  first  three  steps.    The   solution 

y^   of  step  two  and  the  four  solutions  of  step  three  have  a  common  domain  of 
r 

definition.    Thus   y„  is  a  linear  combination  of  the  four  other   solutions. 

F 

Since   y.  and  y.  are  exponentially  large  in  the  common  domain  y„   must   be   a 
14  p 

linear   combination   of  y.  and  y, .    When  one  matches  the  amplitudes  of   y_, 

2       3  f 

y  ,  and  y, ,  one  sees  71?  ~  y^  "•■  ^^^i'  "^^  '^®  proof  is  done.  Thus,  we  prove 
only  steps  one,  two,  and  three. 

We  start  the  proof  with  the  construction  of  a  comparison  equation.    We 
define 

M^  =  (y^  yI  73  y^).  (All) 

and  we  introduce 

B^  =  (t^  M^)  (M^)'^  (A12) 

ds 
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Jince  det  M   is  limited  by  (A8)  we  know  that  B   exists.   Farther 

-p  M^  =  B^  M^.  (A13) 

as 


so  that  kT  is  a 


fandamental  matrix  of  solutions  for  the  equation 


y'  =  B^  y.  (A14) 


Further,  by  construction  and  the  uniqueness  of  the  asymptotic  expansions  we 
know  that 

P  % 

Ib  -  AI  i  C^/lsl  "  (A15) 

where  both  P  and  Qp  go  to  infinity  together.  We  could  determine  the  rela- 
tion between  P  and  Q  ,  but  it  is  of  no  significance  in  the  proof.  We  may 
rewrite  (Al)  as 

y'  =  B^  y  +  (A  -  B^)y, 


which  can  be  written  as  an  integral  equation 

y  =  y^Cs)  +  /(s)  /*  [M^{s')]"^  [A(s')  -  B(s')]  y(s')  ds',         (A16) 

P    P 

where  y   is  any  solution  of  (A14) ,   that  is  a  linear  combination  of  y  ,   y  , 

73.  and  y^. 

In  order  to  construct  solutions  of  (A16)  we  must  specify  the  lower  lim- 
its of  the  integration  in  (A16),  specify  y-(s),  and  then  demonstrate  that  a 
unique  solution  of  (A16)  exists.  We  would  then  also  like  to  show  that  the 
asymptotic  expansion  of  this  solution  is,  up  to  some  order,  y  (s).  We  will 
only  sketch  the  proof  which  is,  in  every  case,  a  straightforward  application 
of  Picard  iteration.  Picard  iteration  will  also  show  that  the  asymptotic 
expansion   of  y(s)  is,   up  to  some  order  y  (s).    In  order  to  describe   the 
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proof  farther  it  is  convenient  to  introduce  the  notation  for  the  exponential 
factors  in  (A4)-(A7): 


z^   =  exp{is  +  i(b/2)log  s]  =  l/z^  (A17) 

3/2       1/2 
z,  =  exp(cs'"''  +  d  s^''')  =  1/z..  (A18) 

Then,  we  may  identify  the  structure  of  M^(s)  [M^(s')]    as 

Mp(s)  [Mp(s')]~^ 

=  (s')"^^^[z,(s)z  (s')  N^,(s.s')  +  z  (s)z  (s')  N^,(s.s')l 

-1/4  P  p  (A19) 

+  (ss')   '^[z,(s)z.(s')  n;^(s.s')  +  z.(s)z,(s')  N;-(s.s')]. 
3     4       34  4     3        43 

where   the  four  by  four  matrices  N^ ,  are  all  bounded  in  Im  s  <  0  and  Im   s'< 
'  ij  -  - 

0,  |s|  >.  1,  Is'  I  >.  1,  so  that 

|N^  (s.s')l  i  Lp.  (A20) 

We  now  carry  out  the  second  step  of  the  proof.    We  restrict   ourselves 

to   the   sector  -5n/6  <_  arg  s  ^  0,   and  we  set  y   in  the   integral   equation 

p 
equal  to  y  (s).    We  select  different  end  points  of  integration  for  the  four 

types   of  terms  that  appear  in  the  kernel  as  shown  in  (A19) .    We  expect   to 

p 

show  that  y  ~  y     Thus,   for  the  first  two  types  of  terms  z.(s)  z-(s')  and 

z_($)  z,(s')   we   take  the  path  of  integration  to  be  the  line  Re  s'   =   con- 

p 

stant.  On  this  path,  provided  y  -  y.  ~  z.  the  exponential  factors  take 
their  largest  values  at  the  end  point  s'  =  s,  and  we  may  take  these  exponen- 
tial factors  outside  the  integral,  evaluated  at  s'  =  s.  Now  z  is  exponen- 
tially large  in  the  sector  -n/3  <  arg  s  <.  0  and  exponentially  small  in  the 
sector  -n  <  arg  s  <  -n/3.  Where  z,  is  large  z.  is  small,  and  conversely. 
The  integrals  containing   z  (s')  start  at  infinity  in  the  sector  where  z.  is 

small,  and  correspondingly  for  z  (s').   It  is  again  easy  to  adjust  the  paths 

4 

of  integration  such  that  the  exponential  factors  are  largest  at  s'  =   s,   so 
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that  these  factors  can  be  taken  out  of  the  integral.  We  now  define  a  norm 
of  a  function 

II  yll  „  =  sup  II  y(s)  z,  (s)||  .  (A21) 

^   IsllR        ^ 
-5n/6<.arg  s<.0 

It  then  easily  follows  that 

16  L     -Qp+3/2 

l|r(s)  f^    [»r(s')]  ^[A(s')-B*^(s')]  y{s')  ds'll   < R 

(Qp-3/2) 

(A22) 

Thus,   for  any  P  sufficiently  large  that  Qp  >  3/2,   there  is  a  radius  R  such 

that   the  integral  operator  defines  a  contractive  mapping  for  |s|  >  R  in  the 

sector  -5n/6  <_  arg  s  <^  0.   Hence  the  integral  equation  has  a  unique  solution 

in   functions  of  bounded  norm  (A21).    This  solution  satisfies  the   original 

equation   (Al),   and   may  be  extended  uniquely  to  the  full  lower  half   plane 

Im  s   <.  0.    Further,   the  estimate  (A22)  also  shows  that  to   leading   order 

(actually   order  Qp  -  3/2)  the  solution  has  the  same  asymptotic  expansion  as 

p 

y.   It  is  this  solution  we  call  y_,. 
2  r 

We  do  not  give  the  third  step  of  the  proof  as  it  is  a  repetition  of  the 

above   argument   in  the  sector  -n  <.  arg  s  <^  -2jt/3  -  5,   0  <.  5   <.  n/6.    One 

P    P    P        P 

replaces  y   in  (A16)  successively  by  y  ,   y  ,   y. ,   and  y  ,   and  adjusts  the 

paths  of  integration  to  infinity  such  that  the  exponential  factors  in  the 
integrands  are  the  largest  at  the  end  point  s'  =  s.  An  alternative  proof  of 
the  validity  of  the  expansions  in  this  sector,  in  which  no  expansions  change 
from  large  to  small,  or  vice  versa,  is  given  in  Coddington  and  Levinson. 
The  proof  is,  in  any  case,  quite  straightforward,  and  we  omit  it. 
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Fig.l   Contours  C  ,   C  ,   and  C.  in  the  p-plane  which  generate  the  solutions 
o^,  u_,  and  a  ,  respectively,  of  Eq.(27). 
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Fig. 2   Uniform  field  dispersion  relation  for  the  magnotesonic  waves  near  the 
fundamental  ion  resonance. 
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Fig. 


3  Propagation  characteristics  for  waves  across  the  third  ion  cyclotron  resonance 
harmonic  for  T  =  10  keV,  and  for  various  values  of  X^  and  L  =  Lw/c  =  1 ,  as  a 

function  of  n,,.   A^  ,  A^,  A^  denote  the  absorption  coefficients  for  LFI ,  HFI, 

and  Bernstein  mode  incidence,  respectively.   P-^,  and  R^  are  the  reflection 

coefficients  for  LFI,  and  BMI,  R^  and  T^^  denote  the  reflected  and  transmitted 

mode  conversion  coefficients. 
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Fig. 4   Absorption   contours   for  the  third  harmonic  resonance  for  T  =  5   keV 
and  L  =  1  in  the  (Xg-nn  )  plane. 
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Fig.5   Absorption   contoars   for  the  third  harmonic  resonance  for  X  =  1000 

H 

and  L  =  1  in  the  (T,n..  )  plane. 
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Fig. 6   The  effect  of  the  machine  size,   or  L,  on  the  absorption  coefficient. 


for  X   =  1000,  T  =  10  keV,  for  the  third  harmonic  resonance. 
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Fig. 7  Absorption  contours  for  T  =  10  keV,   L  =  .95.   in  the  (Xg.^ii  )   plane 
for  magnetosonic  waves  across  the  fundamental  ion  resonance  layer. 


-62- 


'II 


HFl 

LFI 


T=15keV 
L  =  055 


2000 


4000 


6000 


8000 


1000C 


H 


Fig. 8   Absorption   contours  for  T  =  15  keV,   L  =  .95.   in  the  (Xg.n||  )  plane 
for  magnetosonic  waves  across  the  fundamental  ion  resonance  layer. 
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Fig. 9  Absorption  contours  for  the  ion-ion  hybrid  resonance  in  the  minority 

density  and  parallel  wavenumber  plane  for  a  PLT  example  with  B  =  3T, 

L  =  130  cm,  T_  =  T_  =  2  keV,  n  =  3  XlO^^cm"  . 
u    u  o 
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Fig. 10  The  same  as  Fig. 9  but  for  a  JET  example  with  B  =  3.45T,  L  =  3m,  T^  = 


T„  =   5  keV.  n   =3.5  XlO^^cm  ^, 
H  e 
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Fig. 11  The  same  as  Fig. 9  but  for  a  CIT  example  with  B  =  10  T.   L  =  122   cm, 
Tgg   =  T^  =  10  keV,   %  =  5  XIO   cm  '.        Dashed  lines  denote   level 
lines   for   the  energy  absorption  by  tritium  for  the  low   field   side 
incident  fast  mode. 
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Fig.l2  Power  deposition  profiles  for   He  and  T  in  the  CIT  example  of  Fig. 11, 


for  n   =  0.005  n   and  kn  =  2  1/m. 
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